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ABSTRACT 

We present harmonic decompositions of the velocity fields of 19 galaxies from THINGS (The H I Nearby 
Galaxy Survey) which quantify the magnitude of the non-circular motions in these galaxies and yield obser- 
vational estimates of the elongations of the dark matter halo potentials. Additionally, we present accurate 
dynamical center positions for these galaxies. We show that the positions of the kinematic and photometric 
centers of the large majority of the galaxies in our sample are in good agreement. The median absolute ampli- 
tude of the non-circular motions, averaged over our sample, is 6.7 kms -1 , with ~ 90 percent of the galaxies 
having median non-circular motions of less than ~ 9 kms -1 . As a fraction of the total rotation velocity this 
translates into 4.5 percent on average. The mean elongation of the gravitational potential, after a statistical 
correction for an unknown viewing angle, is 0.017 ±0.020, i.e., consistent with a round potential. Our derived 
non-circular motions and elongations are smaller than what is needed to bring Cold Dark Matter (CDM) sim- 
ulations in agreement with the observations. In particular, the amplitudes of the non-circular motions are not 
high enough to hide the steep central mass-density profiles predicted by CDM simulations. We show that the 
amplitudes of the non-circular motions decrease towards lower luminosities and later Hubble types. 

Subject headings: galaxies: kinematics and dynamics — dark matter — galaxies: fundamental parameters — 
galaxies: ISM — galaxies: spiral — galaxies: dwarf 



1, INTRODUCTION 

Dark Matter is one of the fundam ental ingredients o f our 
current cosmological paradigm (e.g.. ISpergel et al1l2007l) . It 
is the dominant mass component in most galaxies, and as such 
determines the properties and evolution of these objects. Mea- 
suring or modeling the properties of these dark matter halos 
can thus help improve our knowledge of galaxy evolution and 
its relation with cosmology. Dark matter cannot be observed 
directly, so its characteristics must be derived observationally 
by studying the gravitational effect it has on tracer particles 
(be it atoms or photons). Alternatively, its properties can be 
studied using sophisticated computer simulations. The latter 
have now converged on a picture where dark matter consists 
of massive, non-relativistic particles, also called Cold Dark 
Matter (CDM). 

Simulations adopting this framework find galaxy ha- 
los to be tri-axial (e.g., lFrenketal.ll 1988UDubinskil 1 1994t 
lHayashi et alJ20 04. 2007: lCapuzzo-Dolcetta et alJ2007l) . and 
to have a characteristic mass density profile. Its most dis- 
tinghuishing feature is the steep increase in dark matter 

?ed as a power - 
afl [T996l Il997h . 
Much effort has been dedicated to obtaining observational 



density towards the center, usua lly described as a power - 
law p ~ r a with a < — 1 (e.g., iNavarro et 
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confirmation of these predictions, concentrating mainly on 
dark matter dominated galaxies. Most observational de- 
terminations of the distribution of dark matter in these 
galaxies (as deduced from their kinematics) seem, how- 
ever, to indicate that it is not characterized by a r" 1 "cusp", 
but by a centra l kpc-sized "core" with a constant dark 
matter d ensity dMcGaugh & de Blokl 119981: Ide Blok et al.1 
1 2001 albl: iMarchesini et all l2002t Ide Blok & Bosmal 1 2002 1 : 
Gentile et all I2004t ISimon et al.ll2003l 120051: Ide Blokl 2005 : 
Zackrisson et al.ll2006t IGentile et . alj|2007t ISpano et alj|2008 : 
iKuzio de Narav et alj|2006l 120081) . 

Whether this core is real or a result of unrecognized sys- 
tematics in the observations has been the subject of some de- 
bate. For example, it is conceivable that the photometric and 
kinemati c centers could be physically offset from each other 
(see, e.g. [Matthews & Gal lagher 2002), which would lead to 
an artifical flattening of the inferred mass de nsity distribu- 
tion if the "wrong" center is chosen (see, e.g., ISwaters et al.l 
l2003al) . A second effect which has been discussed by many 
authors is the effect of non-circular m otions. These are a 
natural consequence of t r i-axial halos (lHayashi et all l2004t 
lHavashi & Navarro! 120061: lHavashi et al]|2007l) . but can also 
be caused by ba rs, spiral arms, or other perturbations t o 
the potential (cf. iRhee et all l2004t IValenzuela etail 120071) . 
As an example of how t hey could affect the observations, 
lHayashi & Navarrol(l2006l) show that the ellipticity of the dark 
matter halo can induce large non-circular motions in the inner 
parts of galaxies (up to ~ 15 percent of the total rotation ve- 
locity), thus making an intrinsic cuspy density profile appear 
cored. Using simulated rotation curves, de BloketaL|(2003) 
show that non-circular motions of the order of ~ 20 kms" 1 
over a large fraction of the disk are needed to make them con- 
sistent with CDM halos (see also IGentile et al.ll200H who 
reach similar conclusions derived from actual observations). 

As mentioned above, most of these studies have concen- 
trated on dark matter dominated galaxies, and more specifi- 
cally on extremely late-type spirals as well as late-type dwarf 
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galaxies (a group of galaxies commonly known as Low Sur- 
face Brightness [LSB] galaxies). The baryonic components of 
these galaxies, being a minor contribution to the total mass, 
can be used in these systems as a tracer for the dark mat- 
ter component. However, a question which naturally arises 
is how representative these galaxies are for the galaxy pop- 
ulation at large. The shallow potentials and low mass sur- 
face densities potentially allow the possibility that their pho- 
tometric and kinematic centers have particularly large offsets, 
and/or that the importance of non-circular motions in their 
disks is anomalously high. Systematic observational studies 
of center positions and non-circular motions in late-type spi- 
ral and dwar f galaxies are, however, still comparatively rare 
(but see, e.g ., Matthe ws & Gallaghed2002tlSimon et al.l2003l 
120051; IGentile et al.ll2005l: iKuzio de Narav et alj|2006fl 2008). 

One can also take a broader approach and study non- 
circular motions and center offsets across the Hubble se- 
quence, with the aim of, firstly, tying the existing studies 
of dark matter-dominated galaxies into a larger, established 
scheme of galaxy properties changing with Hubble type and 
luminosity, and, secondly, providing a baseline for future 
studies. 

For this reason, we have embarked on a major study of a 
large sample spanning a major part of the Hubble sequence 
to measure possible offsets between positions of the dynami- 
cal and photometric centers, as well as the importance of any 
non-circular motions. We will be using the THINGS (The 
HI Nearby Galaxy Survey; [Walter et al.ll2008l see Section 2) 
sample which covers a large range in physical properties, and 
enables us to address aspects such as the importance of dark 
matter as a function of Hubble type, the applicability of dark 
matter models, the shapes of rotation c urves, etc. Many of 
thes e topics a r e cove red in some detail in lde Blok etail £2008) 
and lOh et al.1 (l2008h . The THINGS sample furthermore con- 
tains a number of dark matter dominated dwarf galaxies that 
should provide direct clues on the behaviour of the centers and 
the importance of non-circular motions in extremely late-type 
galaxies and enable a comparison with earlier Hubble types. 

Deviations from purely circular motions in galaxies in gen- 
eral can have different causes, including chaotic non-circular 
motions induced for example by star formation (a s investi- 
gated in detail for two THINGS dwarf galaxies in lOh et alj 
2008), or systematic non-circular motions that relate to the 
potential (e.g., spiral arms, tri-axiality of the halo). Here, we 
will focus on the systematic non-circular motions. One way 
to quantify these is to make a harmonic decomposition of 
the velocity field. An extensive descri ption of the harmoni c 
decomposition technique is given in ISchoenmakersI (119991) . 
Several groups have measured no n-circular motions in disk 
galaxies using this technique (e.g..lSchoenmakers et alJll997t 
IWong et al.1 l2004t IGentile et alTl2005l) . However, none of 
these works report non-circular motio n s as large as those 
proposed by lHavashi & Navarro! d2006l) . ISimon et all (120031) 
used a different technique to measure the non-circular mo- 
tions in NGC 2976. They conclude that although the inner 
300 pc of NGC 2976 contain relatively large non-circular mo- 
tions, the slope o f its (cored) density pro file is not significantly 
affected by that. ISwaters et al.l (2003b) on the contrary argue 
that the non-circular motions in DDO 39 might prevent the 
derivation of a well-constr ained density profile. 

As f or the center offsets. IdeBlok et al.l ( l2001al) and lde Blokl 
(2004) show that significant differences between the mass 
density profiles of cored and cuspy halos become prominent 
interior to a radius of ~ 1 kpc. This radius follows directly 



from the halo models and is fairly independent of (realistic) 
assumptions on halo parameters of either model. Center off- 
sets significantly smaller than a kpc are thus unlikely to have 
a major impact on the shapes of the rotation curves. 

A previous study by [Matthew s & Gallagher (2002) directly 
addresses the issue of center offsets using high-resolution 
rotation curves of a sample of 21 extreme late-type, low- 
luminosity, low-surface brightness spiral galaxies. They com- 
pare the systemic velocities as derived from HI and Ha obser- 
vations and in general find excellent agreement (for two-thirds 
of their sample the velocities agree to better than 10 km s" 1 ), 
leading them to conclude that "most extreme late-type spirals 
seem to have well-defined centers in spite of their often dif- 
fuse stellar disks and shallow central potentials." 

In this paper we will quantify the central positions of the 
THINGS galaxies, as well as quantify the strength of the 
non-circular motions in these galaxies, and investigate their 
dependence on global galaxy properties, such as luminosity. 
We will also check wheth er the mass models o f the THINGS 
galaxi es as presented in Ide Blok et alJ (l2008h and lOh et alJ 
(2008) are affected by these non-circular motions. 

The paper is organized as follows. After a brief descrip- 
tion of the sample in §2, we extensively discuss the various 
ways in which one can define the center of a galaxy in §3, and 
present our best determinations of the positions of the true 
kinematic centers of our sample galaxies. These center posi- 
tions are used as inputs for our harmonic decompositions of 
the observed velocities, as described in §4. We illustrate the 
procedure and show the results for one case-study galaxy in 
detail in §5. Results for the rest of the sample are presented 
in the Appendix. Our results and several quality checks are 
presented and discussed in §6, and the conclusions are sum- 
marized in §7. 

2. SAMPLE AND DATA 

For our analysis, we used THINGS (The H I Nearby Galaxy 
Survey), an H I spectral line survey of 34 nearby disk galaxies 
obtained at the NRAO 6 VLA in B, C, and D arrays. THINGS 
contains a wide range of galaxy types, i.e., high and low sur- 
face brightness galaxies, grand design spirals and dwarf ir- 
regulars, barred and non-barred galaxies, all observed at high 
spatial (~ 10" for the natural weighted data cubes) and ve- 
locity (< 5.2 k ms" 1 ) resolution (f or a detailed description 
of THINGS, see IWalter etail l2008). The great advantage of 
THINGS is that all galaxies have been observed, reduced, and 
analyzed in a homogeneous manner. The large range in phys- 
ical properties enables one to study trends with, e.g., Hubble 
type. Furthermore, as THINGS was designed to overlap with 
the SINGS survey (Sp itzer Infrared Nearby Galaxies Survey, 
iKennicutt et al.ll2003l) . we can use its IRAC 3.6 /jm images 
(which give a virtually dust-free view of the stellar disk) to 
directly compare photometric and kinematic centers. 

In this paper, we will study the kinematics of galaxies. 
We therefore use t he sub-sample of THINGS for which 
Ide Blok et al.l (|2008) derived accurate rotation curves. This 
sample includes THINGS galaxies with inclinations larger 
than 40° that are not undergoing obvious strong interactions, 
as well as NGC 6946, which, despite its lower inclination 
proved suitable for the derivation of a rotation curve. Basic 
properties of the galaxies of our sample are given in Table Q] 

6 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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For our analysis, we used the velocity fields which were cre- 
ated by fitting hermite polynomials to the velocity profiles 
from the natural we ighted data cubes (as described in detail in 
Ide Blok eta l. 2008). Additionall y, we made use of the radio 
continuum maps from THINGS (Walter et al. 200 j), as well 
as the 3.6 /im images from SINGS dKennicutt et al.l l2003). 

3. ESTIMATING GALAXY CENTERS 

For a correct appraisal of the non-circular motions in a 
galaxy (as derived from a harmonic decomposition) it is im- 
portant to accurately determine the position of its center (cf. 
ISchoenmakerslfl999h . There are several ways in which one 
can determine the position of the center of a galaxy. As central 
activity in galaxies is likely to coincide with the bottom of the 
galaxy potential well, a central compact radio source is a good 
indicator of the center of a galaxy. Simil arly, nuclear star clus- 
ters c an be used to locate the center (Matthews & Gallagher 
2002). One can also fit ellipses to the surface brightness distri- 
bution as well as construct a tilted-ring model using the kine- 
matic data to obtain additional, independent estimates. 

Ideally, these different determinations should all coincide. 
Large discrepancies can indicate strong disturbances, caused, 
e.g., by spiral arms or strong bars, or a genuine offset be- 
tween the kinematic and photometric center. For the major- 
ity of our galaxies, multiple center estimates can be derived, 
namely kinematic centers derived from the velocity fields, as 
well as photometric centers as deduced from the 20-cm radio 
continuum maps taken from THINGS and from the Spitzer 
3.6 /im images. 

As will be shown in the subsequent analysis, for most 
galaxies the kinematic center position agrees with the posi- 
tion of the radio continuum source and/or the center as derived 
from the IRAC image within the uncertainties. In those cases 
where the positions of the photometric and kinematic centers 
could be determined accurately, we have generally adopted 
one of the photometric centers as our best center position, 
given that the photometric centers usually have smaller uncer- 
tainties. Note that this particular choice affects our determina- 
tion of the non-circular motions. The position of the dynam- 
ical center is to first order that position which minimizes the 
non-circular motions in a tilted-ring fit and adopting this (or 
a variable) center position would thus results in non-circular 
motions slightly smaller than the ones derived here for the 
photometric centers. For most of our galaxies the difference 
will be minimal, but the general effect is important to keep in 
mind. 

Also note that it is still important to evaluate any possible 
differences between the two types of centers, even when both 
are well-defined. This will enable a better understanding of 
possible systematic effects inherent in the method, which is 
important for the interpretation of the results from galaxies 
where the centers may not be as well-defined. 

Below follows a description of the methods used to deter- 
mine the various center estimates. An application of these 
methods is given for one "case-study" galaxy in Section [5] 
Detailed and more technical descriptions are then given for 
all galaxies in Appendix lAl 

3.1. Radio continuum 

A nuclear point source in the radio continuum is usually as- 
sociated with a central compact object, which naturally should 
be at (or very close to) the bottom of the potential well of 
the galaxy. For all galaxies which show such a source in the 
THINGS radio continuum maps, its position was determined 



by fitting a Gaussian to the central source. The uncertainties 
for the center positions estimated in this way are all similar 
(< 1") and we therefore do not show individual uncertainties 
for center estimates deduced from the radio continuum. 

3.2. Spitzer/IRAC 3.6 fitn image 

We m ake use of the high-r esolution 3.6 /im images from 
SINGS dKennicutt et alJ l2003) 7 . These allow an almost dust- 
free view of the predominantly old stellar populations, though 
we note that the 3.6 /xm band can also contain some trace 
emission from hot dust, PAHs and AGB stars. 

For all galaxies in our sample, we determined the central 
3.6 /im position by fitting ellipses using the GIPSY 8 task ELL- 
FIT, taking care that the ellipse fits were not affected by small 
scale structures. For those galaxies which also show a well- 
defined nuclear source in the 3.6 /im image, we additionally 
derived the central position by fitting a Gaussian to the cen- 
tral source. The two different center estimates generally agree 
very well. However, as the center position derived by fitting 
a Gaussian to the central source is usually better constrained 
as the center from ELLFIT, we only list the former (where 
available) in Table [2] Because of the homogeneous and con- 
sistently small positional uncertainty of less than 1", we do 
not list these here. 

3.3. Kinematic center 

In addition to the photometric centers mentioned above, we 
also derive kinematic centers using the GIPSY task ROTCUR. 
This task fits a set of tilted rings of a given width to the ve- 
locity field of a galaxy and determines their central positions, 
rotation and systemic velocities, inclinations and position an- 
gles. We use the best available center position (i.e., a central 
continuum source, if present, otherwise a nuclear source in the 
3.6 /im image and as a last resort the center as derived using 
ELLFIT) as an initial center estimate for ROTCUR and make 
a fit with all parameters left free (including the center). By 
averaging the central positions over a radial range unaffected 
by spiral arms or other large-scale disturbances, we derive the 
position of the kinematic center for each galaxy. 

The determination of the positions of all centers are de- 
scribed in Appendix [A] and all center estimates are summa- 
rized in Table [2] where our adopted best center positions are 
shown in bold face. 

4. HARMONIC DECOMPOSITION 

We perform a harmonic decomposition of the velocity fields 
by decomposing the velocities found along the tilted-rings 
into multiple terms of sine and cosi ne. 

Following ISchoenmakersI (1 19991) . we describe the line-of- 
sight velocity, vi os , as: 

N 

ViosW = v sys (r) + ^c„,(r) cos mip + s m (r) sin m ip, (1) 

m=l 

where N is the maximum fit order used, r is the radial dis- 
tance from the dynamical center, tp is the azimuthal angle in 
the plane of the disk, and v sys is the th order harmonic com- 
ponent, co. Initial tests showed that a decomposition of the 
velocity fields up to third order (i.e., N = 3) is sufficient to 

7 A small number of galaxies in our sample were not part of SINGS. For 
these galaxies, the data were retrieved from the Spitzer archive. 

8 G IPSY, the Groningen Image Processing SYstem ivan der Hulst et all 

GUI 
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TABLE 1 

Basic properties for the galaxies in our sample. 



Name 


D 


?2S 


M B 


i PA 


Vtot 


Ar 


Mm 


Type 




Mpc 


kpc 


mag 


O O 


kms~' 


It 


10 s M 




(1) 


(2) 


(3) 


(4) 


(5) (6) 


(V) 


(8) 


(9) 


(10) 



NGC 925 


9.2 


14.2 


-20.04 


66 


287 


115 


3.0 


45.8 


7 


NGC 2366 


3.4 


2.2 


-17.17 


64 


40 


55 


6.0 


6.5 


10 


NGC 2403 


3.2 


7.4 


-19.43 


63 


124 


135 


4.0 


25.8 


6 


NGC 2841 


14.1 


14.2 


-21.21 


74 


153 


260 


5.0 


85.8 


3 


NGC 2903 


8.9 


15.2 


-20.93 


65 


204 


190 


7.0 


43.5 


4 


NGC 2976 


3.6 


3.8 


-17.78 


65 


335 


80 


3.5 


1.4 


5 


NGC 3031 


3.6 


11.6 


-20.73 


59 


330 


200 


6.0 


36.4 


2 


NGC 3198 


13.8 


13.0 


-20.75 


72 


215 


150 


6.0 


101.7 


5 


IC 2574 


4.0 


7.5 


-18.11 


53 


56 


70 


6.0 


14.8 


9 


NGC 3521 


10.7 


12.9 


-20.94 


73 


340 


210 


6.0 


80.2 


4 


NGC 3621 


6.6 


9.4 


-20.05 


65 


345 


140 


6.5 


70.7 


7 


NGC 3627 


9.3 


13.9 


-20.74 


62 


173 


133 


5.0 


8.2 


3 


NGC 4736 


4.7 


5.3 


-19.80 


41 


296 


120 


5.0 


4.0 


2 


DDO 154 


4.3 


1.2 


-14.23 


66 


230 


48 


6.5 


3.6 


10 


NGC 4826 


7.5 


11.4 


-20.63 


65 


121 


150 


5.0 


5.5 


2 


NGC 5055 


10.1 


17.4 


-21.12 


59 


102 


190 


5.0 


91.0 


4 


NGC 6946 


5.9 


9.8 


-20.61 


33 


243 


220 


3.0 


41.5 


6 


NGC 7331 


14.7 


19.6 


-21.67 


76 


168 


233 


3.0 


91.3 


3 


NGC 7793 


3.9 


6.0 


-18.79 


50 


290 


130 


6.0 


8.9 


7 



NOTE. — (1): the name of the galaxy; (2): distance as given in Walter et al. (2008); 
(3): radius of the major axis of the galaxy at the fig = 25 mag arc sec' 2 isopho t e leve l, 
taken from LEDA; (4): absolute fi -band magnitude as given in Walter etiD <2008|) ; 
(5): avera ge inclination as give n in[de Bloket al. (2008j); (6): average position angle 
as given in de Blok et al. 1 2008); (7): total rotation velocity (used for the normalization 
of the n on-circular motions ); (8): adopted spacing of the tilted-rings; (9): H I mass as 
listed in Walter et al. I200§j); (10): morphological type from LEDA 



capture most of the non-circular signal, as is described in Sec- 
tion |6AT] 

The usual description of the apparent velocity, under the 
assumption of purely circular motion can be retrieved by only 
including m = and m = 1 terms in Eq.Q] i.e., 

ViosO) = v sys (r) + ci(r) cos ip + ij(r) sin V, (2) 

and by ignoring streaming (radial) motions (i.e., S\ = 0). The 
circular rotation velocity corresponds therefore to c\. Note 
that the dependence on inclination is included in the c m and s m 
terms. For the following discussion it is worthwhile to repeat 
a few rules o f thumb which apply t o harm onic decompositions 
as giv en in ISchoenmakers et all (LL997) and ISchoenmakersI 
(fl999l) : 

(1) A perturbation of the gravitational potential of order m 
will cause m+ 1 and m — 1 harmonics in the velocity field (so 
an m = 2 two-armed spiral component will cause m = 1 and 
m = 3 harmonics in the velocity field). 

(2) Perturbations in the gravitational potential are indepen- 
dent and can therefore be added. The same holds for velocity 
perturbations. 

(3) The elongation of the potential e pot in the plane of the 
disk of the galaxy can be calculated at each radius as follows: 

l+2q 2 + 5q 4 

e pot sin2^ 2 = («3-«i) jz jr~, (3) 

ci(l-^ 4 ) 

where q = cos i. The only remaining unknown quantity is if 2, 
the unknown angle in the plane of the ring between the minor 
axis of the elongated ring and the observer. 

(4) Velocities induced by a global elongation of the poten- 
tial will result in a constant offset in e pot sin 2^2- Velocities 
induced by spiral arms occur on much smaller scales, and 
will therefore only lead to perturbations ("wiggles") around 
this offset. 



(5) If the fitted inclination is close to the intrinsic inclination 
of the disk, then C3 = 0. Small offsets of a few kms" 1 result in 
only small (1-2 degree) inclination offsets. 

We use the GIPSY task RESWRI. This task performs a 
tilted-ring fit assuming circular rotation, creates a model ve- 
locity field, subtracts this from the original velocity field, and 
does a harmonic expansion of the residuals. RESWRI does 
not down-weight velocities along the minor axis as is usually 
done in standard rotation curve analysis, and thus emphasizes 
the non-circular motions. For the width of the annuli (cf. Col 
(8) of Table [TJ, we chose half the beam width; neighboring 
rings are thus not independent. 

We calculate the quadratically added amplitude ("power") 
for each order of the harmonic decomposition using 

Ai(r) = y^W, (4) 

for m = 1 (note that c\ corresponds to the circular velocity and 
is not included in the calculation of the amplitude of A\(r)), 
and 

AJr) = ^c 2 Jr) + s 2 m (r), (5) 

for m > 1 . 

Additionally, we calculate the quadratically added amplitude 
of all (i.e., up to N = 3) non-circular harmonic components 
("total power"): 

A,(r) = yj s\{r) + c 2 (r) + s 2 (r) + c\{r) + s 2 (r). (6) 

The radial variation of A m {r) and A r (r) can be checked for 
coincidence with visible features in the galaxies. To quantify 
the power of the non-circular motions in a compact way, we 
also use A m , the median of A m (r) defined for each value of m, 
and A r , the median of A, (r). We derive two values for A m and 
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A r , one for the entire radial range and one for the inner 1 kpc 
(for those galaxies where sufficiently high signal-to-noise H 
I is present in the inner parts). The 1 kpc value is not arbi- 
trarily chosen but follows from both the core and cusp models 
for reasonable choices of halo parameters. The distinction 
between a cusp a nd a core can be made most clear ly inte- 
rior to this radius dde Blok et alj|2001at Ide Blokll2004l) . and it 
is therefore important to separately quantify the non-circular 
motions in this inner region. Values determined for the entire 
radial range are thus determined over the entire extent of the H 
I disk, except for radii where RES WRI failed to converge prop- 
erly or produce stable results. These radii are usually charac- 
terized by very small filling factors and are mostly found in 
the outer parts of the galaxies. 

In the following we will discuss the absolute non-circular 
motions, as well as the same motions expressed as a fraction 
of the local and total rotation velocity. The latter is defined as 
the rotation velocity of the flat part of the rotation curve, or 
as the maximum rotation velocity if the rotation curve is still 
rising at the outermost point. 

RESWRI creates a residual velocity field which can be used 
as an indicator for those non-circular motions that are not cap- 
tured w ith the harmonic decomposition (as will be shown in 
Section |6.4.11 i. As the residual velocity fields contain only 
values which scatter around a mean zero level, we use abso- 
lute residual velocity fields in our analysis. 

Finally, as mentioned in ISchoenmakersi (1 19991) . the center 
position should be kept fixed during the harmonic decomposi- 
tion. This is because a galaxy which has real, physical C2 and 
52 terms in its velocity field would, in the case of an uncon- 
strained center position, appear to have a center which drifts 
in such a way as to minimize these terms. Moreover, rapidly 
varying center positions at small radii (with offsets larger than 
the relevant ring radii) have no physical basis. 

5. NGC 3 198 — A CASE-STUDY 

Here we present the derivation of the central position and 
the harmonic decomposition in some detail for one galaxy. 
The purpose of this section is to explain the conventions and 
notations used and demonstrate our methods for one galaxy 
in our sample. A complete description and discussion of all 
galaxies from our sample is given in the Appendix. 

5.1. Center estimates 

The 3.6 fjm IRAC image of NGC 3198 shows two well- 
defined spiral arms, emanating from a prominent bulge. The 
central component has a nuclear point source embedded, 
which has a counterpart in the radio continuum. The IRAC 
and continuum centers agree to within 1". For ease of refer- 
ence, and in order to have a compact notation for the cen- 
ter positions, we will in our discussion frequently refer to 
positions with res pect to the pointing center (as listed in 
IWalter et alj 120081) . This pointing center has no physical 
meaning, and merely provides a convenient zero-point, unre- 
lated to any particular choice of the central position. We will 
refer to the offsets from the pointing center as SX (positive in 
the direction of decreasing right ascension) and SY (positive 
in the direction of increasing declination). Both are expressed 
in arc seconds. Note that we list the full coordinates of all cen- 
ter positions we derive in Table [2] (with the position used for 
further analysis shown in bold face). 

In our determination of the kinematic center of NGC 3198, 
we start with an unconstrained ROTCUR fit with the position 



of the radio continuum center (SX= —I'M, SY= -0."1) as an 
initial estimate. Fig. Q] shows the variation of the center posi- 
tion from the ROTCUR fit over the radial range of the galaxy, 
together with our best photometric center. As can be seen 
by the variation of SX and SY, the outer parts (r > 150") of 
NGC 3198 are strongly affected by the spiral arms. For the 
derivation of the dynamical center, we therefore restrict the 
averaging of SX and SY to radii with r < 100" (indicated in 
Fig.[TJby the vertical lines at these radii). The dynamical cen- 
ter derived in such a way is offset from the pointing center by 
5X = -l'M±l'.'l, SY =-0."2±2"l (i.e., to the south-east). 
Averaging SX and SY over all tilted-rings results in a similar 
center position, though with a larger scatter. 

To put our estimates for the center positions in context, we 
show them together with the IRAC and radio continuum map 
(where available) overlaid on the central 150" x 150" of the 
H I total intensity map (cf. Fig. |2j. Also shown are the central 
positions of the individual tilted-rings. As can be seen, all 
center positions agree well within the uncertainties and to 
within one natural-weighted beam (hereafter referred to as 
"the beam") and we therefore adopt the center as derived 
from the radio continuum map as our best center position. 



5.2. Harmonic expansion 

The radial distribution of all fitted parameters from the har- 
monic decomposition of NGC 3198 are shown in Fig. [3] The 
PA of NGC 3 198 rises swiftly within the inner 200", and then 
slowly declines. The inclination varies in the inner parts over 
a range of about five degrees, but shows a steady increase be- 
yond r ~ 450", indicating that the outer disk is warped. The 
C3 term is small for all radii, meaning that the fitted inclination 
is close to the intrinsic inclination of the disk. Although there 
is no global offset from zero for C2 and S2, they show small 
deviations at radii coinciding with the locations of spiral arms 
in the 3.6 /i m image of NGC 3198, or in the total H I map 
presented in lWalter et al.l (120081) . The si and 53 terms are best 
described as wiggles caused by spiral arms on top of a slight 
offset. The amplitudes of all non-circular components (C2, C3, 
*i> *2, S3) are generally only a few kms" 1 . 

In Fig. |U we show derived parameters which were calcu- 
lated according to Eqs.[3]|6] The median amplitudes of the in- 
dividual harmonic components, derived following Eqs.|4]and 
[5] are similar in amplitude, ranging from A m ~ 2 - 3 km s" 1 
(when averaged over the entire radial range), or < 2 per- 
cent of y tot . For the inner 1 kpc, the amplitudes are even 
smaller (A m < 1.5kms _1 ). The distribution of A r (r) (cf. Eq. 
|6]l shows that the amplitude of the non-circular motions is 
A r (r) < 8 kms" 1 for most radii. The median amplitude is 
A r ~ 4.5 kms" 1 when averaged over the entire radial range, 
and A r ~ 1 .5 kms" 1 when averaged over the inner 1 kpc only. 

The elongation of the potential, e pot sin(2<£>2), derived ac- 
cording to Eq. [3] is small over most of the radial range, but 
shows traces of the same spiral arms which cause the variation 
in, e.g., the 53 component. Inwards of r ~ 30", the elongation 
is rather unconstrained — mainly because of the larger un- 
certainty in the derived inclination, which enters into the un- 
certainty in e pot as a fourth power. The weighted mean elon- 
gation of the potential is fairly small with (e pot sin(2</?2)} = 
0.017 ±0.020, and within the uncertainties consistent with 
zero. 

As mentioned in Section [4] RESWRI creates residual ve- 
locity fields, which can be used to quantify the signal which 
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FIG. 1 . — Radial variation of the center positions of the individual tilted-rings derived in an unconstrained fit with ROTCUR. The center positions are given as 
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lines. Our best center position (in this case from the radio continuum) is indicated by the dashed horizontal line, which is however barely distinguishable from 
the kinematic estimate because of the good agreement between these two. 



was not captured in the harmonic expansion. For NGC 3198, 
we find a median absolute value of ~ 2.6 km s" 1 , showing 
that a harmonic decomposition up to third order is capable of 
capturing most non-c ircular motions. This is addressed more 
fully in Section I6AT1 

The results from the harmonic decompositions are summa- 
rized (for all galaxies) in Tableland detailed descriptions for 
all galaxies are given in AppendixlAl 

6. RESULTS 

In the previous section and in the detailed notes on the in- 
dividual galaxies in the Appendix, we show how the center 
positions are derived and present the results of the harmonic 
decompositions. The different center estimates are summa- 
rized in Table [2] where our adopted positions are shown in 
bold face. The results from the harmonic decomposition are 
summarized in Table [3] 

In this section, we will put these results in context and ad- 
dress some of the a strop hysical questions discussed in the In- 
troduction. Section l6"Tl will deal with the results of the center 
estimates, whereas Section 16.21 deals with the results of the 
harmonic decompositions. In Section 16.31 we compare the 
latter with predictions from CDM simulations. In Section |6~4"l 
we present and discuss several consistency checks which we 
applied to test our methods. 

6.1. Quality of the center estimates 

In this section, we discuss the quality and reliability of our 
center positions. An intrinsically cuspy density profile can be 
mistaken for a flat, constant-density one, if possible offsets 
between the kinematic and the photometric center of a galaxy 
are ignored. Using the photometric center to derive a rotation 



curve or a mass model in the presence of such an offset will 
result in a less steep rotation curve and density profile. A 
potential cusp could then appear as a core-like density profile. 

As we have determined the kinematic and photometric cen- 
ters of the galaxies in our sample, we can directly test if such 
offsets exist. This is shown in Fig. [6] where we show the off- 
sets between the dynamical and photometric centers for the 
galaxies in our sample. The offsets are generally less than 
^0.1 kpc, with only two galaxies reaching offsets of ~ 0.3 
kpc. 

These results are simil ar to those derived by 
Matth ews & Gallagherl d2002l) for a sample of 21 ex- 
treme late-type spiral galaxies. They compare the systemic 
velocites as derived from the rotation curves of these galaxies 
with the central velocities of the corresponding global HI 
profiles. For over two-thirds of their sample they find offsets 
less than 10 km s" 1 and from this they conclude that their 
sample galaxies, despite their late-type, diffuse stellar disks 
and shallow potentials, do have well-defined centers. A 
similar comparison of the systemic velocities and central HI 
profile velocities (measured at 20 percent of the peak level) of 
our sa mp le galaxies (using the values as listed in [Walter et alj 
2008 and de Bl ok et ai1 l2008) shows a similarly good agree- 
ment with an absolute velocity difference AV = 1 .7 ± 1 .7 km 
s" 1 , the largest offset being 6.2 km s" 1 . 

We also show the distribution of offsets between the kine- 
matic and the best (photometric) center estimates for the 
15 galaxies in our sample with well-constrained photometric 
centers. For 13 out of the 15 galaxies shown in Fig. [5] the dy- 
namical center differs by less than the size of one beam from 
the best (photometric) center — and for 10 galaxies the agree- 
ment is even better than half the size of the beam. Only two 
galaxies (NGC 3627 and NGC 6946) show moderate offsets 
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FIG. 2. — This figure shows the inner 150" X 150" of NGC 3198. The integrated natural-weighted H I map is shown in grayscales. The beam is indicated 
in the bottom-left corner. The black contours show iso-velocity contours from the hermite velocity field. The thick black contour is shown at 663 kms~' and 
represents the systemic velocity given in de Blok et al. 1 2008). The other velocity contours are spaced by 25 kms~' . The thin black contours overlaid on the thick 
white contours represent the 3.6 /im IRAC image. They are drawn at 2, 5, 10, 20, and 50 percent of the maximum intensity level. The white contours represent 
the THINGS radio continuum map, and are drawn at 10, 20, and 50 percent of the peak intensity. The black dots indicate the individual center positions from 
ROTCUR and the black cross represents the derived dynamical center together with its uncertainty. The derived center from the 3.6 fim image is shown as a gray, 
filled triangle, whereas the one from the radio continuum is shown as a black, open triangle. Inset: To better highlight the different center estimates, we show an 
inset of the innermost 18" in the upper-right corner. For clarity reasons we do not show the H I grayscale and velocity field contours here. The contours from 
the 3.6 fim image are shown in black and are given at the same intensity levels as in the main plot. The same holds for the radio continuum contours, which are 
shown here in gray. The individual center estimates from ROTCUR are shown as small crosses. In the inset, the beam is indicated by the thick black dashed ellipse 
and is centered on our best center position. 



of between one and two beam sizes. Given that the former 
galaxy is extremely asymmetric (due to its interactions in the 
Leo group) and the latter has a low inclination which makes 
fitting tilted-rings more difficult, these outliers can be under- 
stood. 

The strength of the agreement we find does depend on how 
well-constrained the dynamical center estimates are. A tightly 
constrained kinematic center position which agrees with a 
well-determined photometric center allows one to draw strong 
conclusions. A weakly constrained kinematic center less so, 
as the increased uncertainties allow agreement with a whole 
range of photometric centers, as long as they are located 
somewhere near the kinematic center. To test the strength of 
our conclusions we therefore determine for each galaxy the 
offsets between the centers of individual tilted-rings and the 
best center estimate. Because of the sparsely filled rings in the 
very outer parts, as well as the presence of warps and asym- 



metries, we restrict this analysis to the data points in the inner 
half of each galaxy (i.e., those with r < 0.5 r max ), with r max as 
listed in Table [3] In Fig. [7] we show a histogram of the dis- 
tribution of these offsets in terms of beam size for all galaxies 
in our sample. Approximately 50 percent of the center posi- 
tions of the ^1000 individual tilted-rings differ less than one 
beam from our best center position. Another 25 percent show 
an offset between one and two beam sizes and only a small 
fraction shows large offsets. Note that the center positions of 
some of the tilted-rings can be affected by the presence of spi- 
ral arms and other features. Large offsets for individual rings 
do therefore not necessarily imply intrinsic offsets between 
center positions, but need to be regarded within the context of 
the results for the whole galaxy. The radial variation of the 
center positions are shown in Fig.Q]for NGC 3198 and in the 
top panels of Figs. I25l42l for the other galaxies in our sample. 
Our results show that for the large majority of the galaxies 
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studied here, the photometric and kinematic center positions 
agree wit hin th eir uncertainties and also to within a beam. See 
also Sec. l6.4.4l where we evaluate the impact of a shift in cen- 
ter position on our results. 

6.2. Results of the harmonic decomposition 

Here, we discuss the results of the harmonic decomposi- 
tions in a more general way. For a description of the results for 
individual galaxies we refer to Appendix lAl Because of th e 
low filling-factor of its velocity field (cf. lde Blok et al.ll2008b . 
we did not perform a harmonic decomposition for NGC 4826. 

In Fig. |8] we show plot the harmonic decomposition re- 
sults as a function of absolute magnitude, Hubble type and 



dark matter dominance. In all panels, we distinguish be- 
tween barred and non-barred galaxies. The Hubble types are 
listed in TableQ]and the bar classifications are based on NED. 
The dark matter dominance is expressed as the ratio of total 
baryonic mass to total dynamical mass M\, m /M tot . The bary- 
onic mass is c alcula ted using the mass models presented in 
Ide Blok et al.l (120081) . where we have used the Kroupa IMF 
results (their Table 4). The total dynamical mass is calculated 
using the radius of the outermost point of the rotation curve 
and the velocity is taken from Table Q] We note that for NGC 
3627 we cannot show M\, m /M tot as the non-circular motions in 
this galaxies were too large for lde Blok et al.l (120081) to derive 
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TABLE 2 

Center positions for the galaxies in our sample. 



Name kinematic (ROTCUR) center 3.6 /im center radio continuum center 
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NOTE. — The uncertainties in the ROTCUR centers are given in units of seconds (for right ascension) and 
arcseconds (for declination). The center position chosen for subsequent analysis is shown in bold face. 

"Derived by fitting ellipses with ELLFIT 

*The cente r coordinates of IC2574 were derived using a bulk velocity field cleared of non-circular motions as 
presented in lOh et alj 120081) 
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FIG. 5. — The histogram shows the offset between the dynamical and the 
best (photometric) center for those 15 galaxies in our sample which have a 
well-constrained photometric center position. The offset is shown in terms 
of beam size. For 13 of the 15 galaxies, the different center estimates agree 
to within one beam (typical beam size: ~ 10"), showing that there is no in- 
dication for a genuine and general offset between kinematic and photometric 
centers in our sample. 



a mass model. 

In the upper panels of Fig. [8] we show the results forA r . It is 
clear that the magnitude of the non-circular motions decreases 
towards lower luminosities and later Hubble types (top panel 
of Fig. A more tentative trend is visible for M\,. M /M tot 
where the dark-matter dominated galaxies have on average 
the lowest non-circular motions. When these motions are ex- 
pressed as a fraction of V tot , we see that their contribution is 
roughly independent of luminosity, Hubble type or dark mat- 
ter dominance (Fig. [8] middle panel). For 16 out of the 18 
galaxies, the non-circular contribution is smaller than 6 per- 
cent of y tot and it has a mean value (averaged over the en- 
tire sample) of 4.5 percent. Both in the top and in the mid- 
dle panel of Fig. [8] barred galaxies do not stand out, except 
for NGC 3627. This is different if we look only at the non- 
circular motions in the inner 1 kpc (Fig. [8] bottom panel). As 
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position of the photometric center as our best center estimate. Open circles 
represent galaxies for which we adopted the position of the kinematic cen- 
ter. Open squares indicate the two galaxies for which no reliable photometric 
center position could be derived. 



expected, most barred galaxies have rather high non-circular 
motions in the innermost region. Note that some galaxies are 
not shown in these diagrams, as we have no data above a de- 
tection limit of 3 a in their inner 1 kpc. 

As a comparison, we also calculated the amplitude of the 
non-circular motions within the inner 2, 3, 4, and 5 kpc, using 
two different methods. The approach shown in Fig. [9] mea- 
sures the non-circular motions by taking the median of the 
A r (r) values within rings of 1 kpc width, as measured over 
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TABLE 3 

Derived quantities from the harmonic decompositions. 
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NOTE. — (1): the name of the galaxy; (2): A,-, the median of the quadratically added 
amplitude of the non-circular motions, averaged over the entire radial range. Uncertainties 
indicate the lower and upper quartile; (3): same as (2) but averaged over the inner 1 kpc 
only; (4): the percentage the non-circular motions contribute to the total rotation velocity; 
(5): the weighted mean elongation of the potential; (6): the median of the absolute residual 
velocity field after the harmonic decomposition; (7): maximum radius for the averaging of 
A- and e pc >t. The bottom two rows contain the mean values and their rms over the entire 
sample. For e po t, this value represents the weighted mean. 

"This value has no upper or lower quartile as the inner 1 kpc of NGC 3198 contain only 
one data point. 

dial range over which the amplitudes are averaged (0 < r < 1 
kpc, < r < 2 kpc, < r < 5 kpc). In both cases, we 
have divided our sample into three absolute magnitude bins. 
It is apparent from the two figures that galaxies with low lu- 
minosity have the lowest amplitudes of non-circular motions 
in the inner parts, regardless of the chosen method. Since the 
approach shown in Fig [9] measures the non-circular motions 
more locally than the method shown in Fig[l0] it is also more 
affected by local features like effects of star formation. For in- 
stance, the two relatively high amplitudes in the 4 < r < 5 kpc 
bin of the bottom panel of Fig. |9]belong to IC 2574 and NGC 
2366. The former has a supergiant shell at those radii (see 
Walter 'etail 11998b . whereas th e latter has a large star form- 
ing region in its outer parts (see lOh et al.ll2008l for an in-depth 
analysis of the small-scale non-circular motions in these two 
galaxies). The only other low-luminosity galaxy in our sam- 
ple which extends out to 5 kpc radius (DDO 154) is com- 
pletely quiescent in contrast. In the analysis shown in Fig.flOl 
these effects have "averaged out", due to the larger area used 
for the averaging. 

However, no matter how they are binned, in absolute terms 
the measured amplitudes of the non-circular motions are 
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FIG. 7. — The histogram shows (for all galaxies in our sample) the off- 
set between the individual tilted-ring center estimates and the corresponding 
best center position in terms of beam size. The typical beam size is ~ 10" 
(this corresponds to a physical size of 170 pc — 750 pc). Only tilted-rings 
inwards of r = 0.5 r max were included. The histogram demonstrates that our 
dynamical center estimates are well-constrained and in good agreemend with 
the photometric center positions. 



radii < r < 1 kpc, 1 < r < 2 kpc, 4 < r < 5 kpc. Our 
second approach, shown in Fig. [10] simply increases the ra- 
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FIG. 8. — The non-circular motions measured in different ways for our 
entire sample. Barred galaxies are shown as rilled circles, non-barred as 
open circles. Top panel: Median amplitude of the non-circular motions (A r ) 
on a logarithmic scale vs. absolute magnitude (left), Hubble type (center) and 
baryonic to total mass ratio (right). Error bars indicate the upper and lower 
quartile of the radial distribution of A r . Middle panel: The percentage the 
non-circular motions contribute to the total rotation velocity. Bottom panel: 
Like top panel, but the averaging was restricted to the inner 1 kpc of the 
galaxies. Note that for five galaxies, the data do not show significant amounts 
of H I in the inner kpc, and they are therefore not shown in the two bottom 
panels. 




20 




30 



20 




M n > -18.5 




2 3 4 5 

mean radius (kpc) 

FIG. 9. — The amplitude of the non-circular motions within rings of 1 kpc 
width (i.e., < r < 1 kpc, 1 < r < 2 kpc, 4 < r < 5 kpc) for each galaxy 
(if data available). All measurements of a specific galaxy are connected 
by a line. Filled symbols correspond to barred galaxies, open symbols to 
non-barred galaxies. We divide the plot into three subpanels. The top 
panel contains all galaxies with Mb < —20.5, the middle panel those with 
-20.5 < M s < -18.5 and the bottom panel all galaxies with M B > -18.5. 



small. In order to account for the different rotation veloci- 
ties, we have normalized the amplitudes shown in Figs.|9]and 
[TO] by the local rotation velocity. The results are shown in 
Figs. [TTJ and [12] respectively. For the large majority of the 
galaxies in our sample, the non-circular motions in the inner 
few kpc contribute approximately ten percent to the local ro- 
tation velocity. Three galaxies (NGC 925, NGC 2903, and IC 
2574) contain non-circular motions larger than 20 percent of 
the local rotation velocity. For NGC 925 and NGC 2903, it is 
likely that the large non-circular motions are associated with 
their stellar bars, especially as Figs. [TT1 and [121 indicate that 
large non-circular motions in the central parts are predomi- 
nantly found in barred galaxies. For IC 2574, the large ratio 
between non-circular and circular motions is mainly caused 
by the slow rise of its rotation curve in the inner kpc. Over 
the same radial range the absolute non-circular motions are 
generally less than 2.3 kms" 1 . 

We will now discuss the elongation of the potential (see 
Fig. [T3V The weighted average elongation and its standard 
deviation are (e pot sin(2</32)) = 0.01 1 ±0.013. This elongation 
and the individual elongation measurements still contain an 
unknown viewing angle tp 2 , for which one can apply a statis- 
tical correction by dividing the elongation by the expectation 
value of ip2, 2 /it (middle panel of Fig.[T3l). The average, cor- 
rected elongation is (e pot ) = 0.017±0.020, i.e., consistent with 
a round potential. 

The elongations show neither a trend with absolute mag- 
nitude, nor with Hubble type, nor with dark matter domi- 
nance. If the influence of the disk decreases with baryonic 
mass, and if the dark matter halos are indeed tri-axial (or more 
precisely, ha ve an elliptical potential di stortion in the plane of 
their disk, cf. Havashi & Navarro 2006), one would expect the 
measured elongations to increase with decreasing (baryonic) 
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FIG. 10. — Like Fig. [5] but the amplitude of the non-circular motions were 
averaged within rings of increasing radius (i.e., < r < 1 kpc, < r < 2 kpc, 
< r < 5 kpc). As the region over which the non-circular motions are 
averaged increases outwards, the amplitudes shown are less affected by local 
features such as the effects of star formation. 



mass. However, no such trend is seen in our data. 

One explanation could be that the baryons dominate the to- 
tal potential sufficiently so that the elongation values we mea- 
sure are really the elongations of the potentials of the baryonic 
disks. As the self-interaction (dissipation) of the (gas compo- 
nent) baryons tends to circularize the orbits, any (kinemati- 
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FIG. 1 1 . — Like Fig. [9] but the amplitudes of the non-circular motions 
are normalized by the local rotation velocity. Except for a few galaxies, the 
non-circular motions in the inner few kpc have an amplitude of ~ 10 percent 
of the local rotation velocity, irrespective of the galaxies' luminosity. Only 
three galaxies have non-circular motions in their centers which are larger than 
20 percent of the local rotation velocity. Note also that large non-circular 
motions (in relative terms) are predominantly found in barred galaxies. 
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FIG. 12. — Like Fig. 1 101 but the amplitudes of the non-circul ar m otions 
are normalized by the local rotation velocity. Analogous to Fig. [10] the in- 
creasing region for the averaging of the non-circular motions causes the am- 
plitudes to be somewhat less affected by local features such as the effects of 
star formation. 

cal) traces of the tri-axial dark matter potential would thus 
be wiped out. With a variation in total M\,. dI /M tot of close to 
an order of magnitude, this means the circularization process 
ought to be extremely effective over a large range of mass ra- 
tios. The relevant quantity here, however, is not the total dark 
matter dominance but the local one. We therefore show in 
Fig. [14] a measure for the local ratio of Mb ar /M tot as a func- 
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FIG. 13. — Top panel: Elongation of the potential vs. absolute mag- 
nitude (left), Hubble type (center) and baryonic to total mass ratio (right). 
The horizontal lines represent the weighted mean and its standard deviation 
((e po t sin(2</32)} = 0.011 ±0.013). Barred galaxies are shown as filled cir- 
cles, non-barred ones as open circles. The distribution of (e pol sin(2</?2)) 
shows neither a trend with luminosity, nor with Hubble type. Middle panel: 
Like top panel, but a statistical correction for the unknown viewing angle 
was applied to the individual elongation measurements and to their mean 
an d standard deviation . The hatched area indicates the CDM predictions 
by Hayashi & Navarro (2006). The elongations of the large majority of the 
galaxies shown here are systematically lower (rounder) than what is predicted 
by CDM simulations, although the measurements for some galaxies have 
large enough error bars to make them marginally consistent with the lower 
end of the CDM predictions. Bottom panel: Like middle panel, but we av- 
eraged the elongation of the potential not over the entire radial range, but 
only out to r < 0.5 r max (^max is tne maximum radius to which our analysis 
extends). This has generally only little effect on the individual elongation 
measurements. The mean value has not changed ((e pot ) = 0.017 ± 0.020 vs. 
(fpot) = 0.016 ± 0.015) and we therefore do not see an indication that the 
elongation of the potentials increase inwards. 
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FIG. 14. — The radial distribution of M hai /M ioi in t he THINGS galaxies , 
derived from the Kroupa IMF mass models presented in de Blok et alj 120081) . 
Light-gray full curves represent galaxies brighter than Mb = 20.5. Dashed, 
dark-gray curves show galaxies with —20.5 < Mb < —18.5. Black curves 
show galaxies fainter than Ms = -18.5. 



tion of r adius, based on the Kroupa IMF mass models pre- 
sented in lde Blok et all d2008[) . We see that whilst for the high 
luminosity galaxies the local importance of dark matter and 
baryons in the inner parts is roughly equal, it drops to values 
of Mbar/Mot ~ 1/3 in the outer parts. For the low luminos- 
ity galaxies the importance of the baryons is, however, close 
to an order of magnitude lower. Again, this implies that the 
circularization process should be equally effective over a large 
range in local mass ratios, and already manage to wipe out the 
tri-axial signal when the baryons contribute only ~ 10 percent 
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to the local mass. Note that in the above we have assumed that 
both stars and gas contribute to the circularization process. It 
is however likely that only the gas component is relevant, as 
the stars will not necessary dissipate as strongly. If this is the 
case, the relevant quantity is the M gas /M tot ratio, which can 
be another order of magnitude lower than the M\, m -/M tot ratios 
discussed here. In summary, attributing the circular poten- 
tials to dissipation requires that the circularization process is 
already effective when the baryons contribute only ^10 per- 
cent to the total mass. Alternatively, the dissipation may not 
be relevant at such low mass ratios, in which case the poten- 
tials are truly round. The well-known dark matter dominance 
of LSB galaxies may imply that similar conclusions also hold 
for them; direct and equally detailed measurements should 
provide the additional information needed to draw definitive 
conclusions. 

6.3. Comparison with predictions from simulations 

One of the clearest ou tcomes of s i mulat i ons is that CDM 
halos are tria xial (e . g.. iFrenket alJ Il988t Dubinski 1 19941: 
Havashi et ail l2004t iMoore et all 12004 iKasun & Evrardl 
20051: lHavashi & Navarro! 120061: lHavashi et alJ 120071) . We 
compare our results with predicti ons based o n cos - 
mological ACDM simu lations. lHavashi etall (120041) : 
lHavashi & Navarrol d2006l) present models where they use an 
elliptical distortion of the gravitational potential (mimicking 
a tri-axial halo) to explain the observed solid-body rotation 
curves of dwarfs and LSB galaxies within a CDM context. 
This distortion affects the orbits and thus also the rotation 
curve. In this picture, the major axis rotation curve depends 
on the viewing angle of the elliptical distortion. This viewing 
angle can be chosen in such a way that the major axis rotation 
curve in the elliptical potential looks identical to what it would 
have been, had the gas been on circular orbits (i.e., with- 
out elliptical distortion and without non-circular motions). It 
would therefore hide the presence of the disturbance due to 
tri-axiality and look like a solid-body rotation curve. 

However, even if the angle is adjusted in such a way, the mi- 
nor axis rotation curve of the tri-axial case would still show 
large streaming motions. These can be up to 15 percent of 
the total rotation velocity. As we show in Fig. [8] the ob- 
served non-circular motions in our sample are generally much 
smaller and never reach 15 percent of the total rotation (see 
also the minor-axis position-velocity diagrams presented in 
Ide Blok etaTJ 120081) . In Figs. [TT] and [121 we normalize the 
amplitudes of the non-circular motions by the local (c i rcular ) 
rotation velocity. At a radius of 1 kpc. lHavashi et al.l (2004) 
report non-circular motions as high as 50 percent of the local 
rotation velocity. At smaller radii, the contribution of the non- 
circular motions is even larger. The results shown in Figs. [TT] 
and[12] however, show that for the large majority of the galax- 
ies in our sample, the non-circular motions are only as high as 
ten percent of the local rotation velocity. Only in one barred 
galaxy (NGC 925) do the non-circular motions reach an am- 
plitude close to 50 percent of the local rotation velocity in the 
inner 1 kpc. The dwarf galaxies (bottom panels of Figs. [TT1 
and[T2b are of most interest here, as these galaxies are the ones 
most dominated by dark matter. The non-circular motions in 
three out of the four dwarfs in our sample (NGC 2366, NGC 
2976, and DDO 154) contribute only ten percent to the local 
rotation velocity. Therefore, the non-circular motions in these 
g alaxies do not signific antl y affect the mass models presented 
in Ide Bloketail d2008l) and lOh et all d2008l) and preclude the 
possibility that their cored density profiles (cf. aforemen- 



tioned papers) are intrinsic cuspy profiles which have been 
artificially flattened by large non-circular motions. The non- 
circular motions in the inner kpc of IC 2574 amount to ~ 26 
percent of the l ocal rotation v e locity , a factor tw o smaller 
than p redicted in lHavashiet all (120041) . We refer to lOhetal] 
(2008) for a full treatment of the IC 2574 non-circular mo- 
tions. 

The gravitational po tential in the simulations of 
lHavashi & Navarrol d2006l) is elongated in the inner parts as 
(e pot ) = 0.2 ± 0. 1 (their Fig. 3). All but two of the individual 
elongation measurements presented here are systematically 
lower than those predictions and consistent with a round 
potential, although some have large enough error bars to 
make them marginally consistent with the lower end of the 
CDM predictions. 

If DM halos do indeed have an elliptical distortion in their 
inner parts, we would expect to find higher elongations there. 
To test this idea, we have also determined the elongations of 
the potential by using only data out to r = 0.5r max (bottom 
panel of Fig. [T3l . The individual elongation measurements 
do not change significantly if averaged only over the inner 
half of the disk. Their weighted mean ( (e pot ) = 0.016±0.015) 
is not significantly different from the value obtained by av- 
eraging over the entire disk ((e po t) = 0.017 ±0.020) and we 
therefore see no indication that the gravitational potential is 
more elongated towards the center of the galaxies. 

One caveat relating to the above discus sion is that the model 
presented in Havashi & Navarro (2006) consists of only a 
dark matter component, while real galaxies contain baryons 
as well. Baryonic processes, such as dissipation in the gas 
component might affect the dynamics within the visible disk 
of observed galaxies, thus invalidating any comparisons. As 
we showed in Sect. 16.21 however, in some of the galaxies in 
our sample the local baryon to total mass ratio is only ~ 10 
percent. If dissipation or circularization of the orbits affects 
the results found for our galaxies, this implies that the 10 per- 
cent of baryons are able to wipe out the dynamical signature 
of the 90 percent of dark matter. If dissipation is really this 
efficient, it would also imply that analyses of, e.g., the scat- 
ter in the Tully-Fisher rel ation in order to put lim its on the 
tri-axiality of haloes (e.g.. lFranx & de Zeeuw|[l992r) are only 
relevant for the elongations of the baryonic disks, but not for 
those of dark matter haloes. In summary, unless the effects 
of baryonic dissipation are extremely efficient, it is difficult 
to reconcile the tri-axialities found in simulations with those 
inferred for the low-luminosity galaxies in our sample. More 
extensive observations and more detailed simulations ought 
to show whether these results hold for the general disk galaxy 
population. 

6.4. Consistency Checks 

In this section, we compare the results of our harmonic de- 
c omposition with the traditional tilted-ring analysis presented 
in lde Blok et alJd2008l) . Additionally, we present quality con- 
trols which show the limitations of our methods, test and val- 
idate our current results as well as indicate room for improve- 
ment. 

6.4.1. ROTCURvi. RESWRI: residual velocity fields 

Here, we compa re the residua l velo city fields from the 
ROTCUR analysis by de B lok et alJ d2008l) with those from our 
analysis with RESWRI. As mentioned in Section[4] the resid- 
ual velocity fields are derived by subtracting a model using 
the final parameter estimates from the original data. For both 
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FIG. 15. — Media n of the abs olute residual velocity field from the rotation 
curve analysis of de Blok et al. 12008) vs. the one from the harmonic decom- 
position presented in this paper. The error bars indicate the lower and upper 
quartile, the solid line represents a one-to-one relation, and the dashed line 
the unweighted least square fit through the data points. The residuals in the 
traditional ROTCUR analysis are all larger than those from the RES WRI fits, as 
expected. The difference between the solid line and dashed line is due to the 
non-circular motions quantified in this paper. 
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FIG. 16. — Median of the absolute residual velocity field from the rota- 
tion curve analysis in de Blok et al. 1 2008) v.v. A,-, the median of the (radially 
averaged) quadratically added amplitudes of the non-circular motions quan- 
tified in this paper. The error bars indicate the lower and upper quartile. 
Galaxies with large residuals in the ROTCUR fits show generally also larger 
non-circular motions. The outlier in the upper part of the panel is NGC 3627, 
which shows large non-circular motions, but only moderate residuals in the 
ROTCUR analysis by de Blok et^al] 120081) . 



types of residual fields, we calculate the median and the lower 
and upper quartile of the absolute residual velocity fields. 
These are compared in Fig. Q3] It is clear that all galaxies 
have larger residuals in the ROTCUR analysis than they have in 
the RESWRI analysis. This is of course to be expected, given 
that ROTCUR considers only circular motion, while RESWRI 
also takes non-circular motions into account. Nevertheless, 
it is a quality test which shows that our results are behaving 
in the expected manner. The largest median amplitude of the 
residual RESWRI velocity fields is 4.5 kms" 1 (for NGC 3521). 
The average value of the sample is 2.9 kms" 1 , which clearly 
shows that a harmonic expansion up to third order is capable 
of capturing the majority of the non-circular motions in most 
galaxies. 

Looking at Fig. [15] it becomes also clear that galaxies with 
small ROTCUR residuals also have small RESWRI residuals. 
Small ROTCUR residuals, however, indicate that a model as- 
suming circular rotation already provides a good approxima- 
tion. Therefore, it is to be expected that these galaxies show 
only small non-circular motions. This expectation is tested 
in Fig. Q2] where we plot the median values of th e abso- 
lute residual velocity fields from lde Blok et al] (120081) vs. the 
median amplitudes A r of the non-circular motions from our 
harmonic decompositions. Most galaxies are located near 
the line of unity, showing that our expectation is indeed cor- 
rect. The only real outlier in Fig. [16] is NGC 3627, which is 
clearly offset from all other galaxies as it has — despite its 
large non-circular motions (A r ~ 28.5 kms -1 ) — only moder- 
ate ROTCUR residuals. The results for the other galaxies, how- 
ever, show that even a simple model considering only circular 
rotation can fit most galaxies quite well. These small resid- 
uals therefore demonstrate (independently from our RESWRI 
analysis) that non-circular motions are generally small. 



6.4.2. Making prior assumptions during the RESWRI runs 

The rotation curves of all gala xies in our sample w ere de- 
rived and discussed in detail by Ide Blok et al.l J2008). The 
general procedure for the derivation of a rotation curve in- 
cludes keeping some parameters fixed in the tilted-ring fits, 
thus reducing the number of free parameters for each indi- 
vidual fit with ROTCUR. In the current analysis we derive all 
tilted-ring parameters in a single fit with all parameters left 
free and do not attempt to correct for the motions induced by, 
e.g., star formation or spiral arms. Thus, it is expected that the 
tilted-ring parameter s derived by us differ to some extent from 
the ones derived by Ide Blok et al.l (120081) . In order to check 
this and to estimate the impact on our results, we compare 
the derived values for our unconstrained harmonic decompo- 
sitions with those from constraine d decompositions w ith PA 
and i fixed to the final values from de B lok et all (120081) . We 
apply this test to two representative galaxies, NGC 3198 and 
DDO 154, a spiral galaxy and a dwarf irregular. 

In the case of NGC 3198 (Fig. [T71 ), most values agree re- 
markably we ll between the c onstra ined and the unconstrained 
fit. Because de Blok et al.l (12008) have assumed a constant 
inclination in the inner parts (in order to compensate for ef- 
fects caused by spiral arms), our unconstrained values differ 
slightly in this region. Note though that the effect on the ro- 
tation curve is negligible. Fixing the inclination has the pre- 
dictable effect that the C3 term reaches values of a few kms" 1 
as this term tries to compensate for the effect the spiral arms 
have on the velocity field. The impact of this, however, is mi- 
nor, as one can see in the distribution of A,-(r) and A m as well 
as in the negligible difference in the derived rotation veloc- 
ity. The elongation of the potential is also largely unaffected 
by constraining i and PA. The only noticeable differences are 
the different-sized error bars. As the inclination contributes 
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to the fourth power to the error of the elongation of the poten- 
tial, the error bars are naturally smaller in the case of a fixed 
inclination. 

For DDO 154 (Fig. [TBI , the differences between the har- 
monic components of the two decompositions are generally 
of the order of ~ 1 kms" 1 . The inclinations in the uncon- 
strained fit differ inwards of r w 80" from the ones in the con- 
strained fit. This is because of the close-to solid-body rotation 
of DDO 154, which makes the simultaneous determination of 
the correct inclination difficult. The impact on the results, 
even in the inner 80", is nevertheless small. For instance, the 
derived rotation velocities are indistinguishable between the 
two decompositions. This shows that the magnitude of the 
non-circular motions is not very sensitive to small variations 
in PA or inclination. 

6.4.3. Hermite vs. Intensity weighted mean VF 

Hermite velocity fields differ from intensity weighted mean 
(IWM) velocity fields when dealing with asym metric (i.e., 
non-G aussian) profiles. It was noted, e.g., by Ide Blok et alj 
(2008), that hermite velocity fields better reproduce the ve- 
locity of the peak intensity of the profile, as their ana- 
lytical function includes an (skewness) term (see, e.g., 
Ivan der Marel & Franxll 1993b . For a detailed overview of the 
different type s of ve locity fields we refer to the discussion in 
Ide Bloket al.l(l2008l 

As we are using the hermite velocity fields, an investiga- 
tion of how using IWM velocity fields affects our analysis, if 
at all, is warranted. We perform unconstrained harmonic de- 
compositions with RESWRI using identical initial conditions 
on the hermite and on the IWM velocity fields of NGC 3198 
and DDO 154. The results are summarized in Figs. [19] and 
l20l respectively. As small differences between the two de- 
compositions are visible, it is clear that how a velocity field is 
created does play a role in the analysis. A close inspection of 
the radial variation of the various quantities shown in both fig- 
ures indicates that the radial variations are more pronounced 
in the case of the hermite velocity field. The choice of the 
velocity field construction method has, however, little impact 
on the derived quantities (see, e.g., bottom panels of Figs. [19] 
and[20"li, and therefore does not change our conclusions in any 
significant way. 

6.4.4. Decomposition under the assumption of an incorrect center 

In this section, we test the influence of an incorrect cen- 
ter position on our results. As test candidates, we have cho- 
sen one galaxy where the position of the center is very well- 
defined (NGC 2841), and one where the center is less ob- 
vious (NGC 2366). For both galaxies, we deliberately shift 
the center positions used in RESWRI by 2", 4", 6", and 10" 
along their major and minor axes. The results of the most ex- 
treme cases (10" offsets, i.e., approximately the beam size) 
are compared with the results from our best center positions 
in Figs. EM] 

For NGC 2366, the effect of an incorrectly chosen center is 
small, irrespective of whether the center is shifted along the 
major (Fig. |2TJ or along the minor (Fig. l22l axis. For off- 
sets along the major (minor) axis, the median amplitude of 
the non-circular motions increases from A,. ~ 3 kms" 1 for our 
best center position to A r ~ 3. 5 kms" 1 (A r ^4kms _1 ). Our re- 
sults for the dwarf galaxy NGC 2366 are therefore unaffected 
by small offsets in the center position, indicating that for a 
rotation curve showing close-to solid-body rotation, a mod- 



est offset from the true center position is not crucial for the 
analysis presented here. 

For NGC 2841, the situation is different. A 10" offset 
along the major axis (Fig. l23l already increases the median 
amplitude of the non-circular motions from A r ~ 7 kms" 1 to 
A r ~ 12kms"'. In the case of an offset along the minor axis 
(Fig.l23l. the increment is even larger (A r ~ 25 kms" 1 ). 

Note that the difference in the harmonic components is 
mostly to be found in the m = 2 term. This is to be expected, as 
a galaxy will appear kinematically lopsided if a center offset 
from the dynamical center is chosen for the harmonic decom- 
position. Also note that, as expected, it is the cosine coeffi- 
cient C2 which shows the largest amplitude in the case of an 
offset along the major axis, and the sine coefficient S2 in the 
case of an offset along the minor axis. 

For NGC 2841 — a galaxy with a steep and subsequently 
flat rotation curve — our results are sensitive to the chosen 
center position and an offset center would clearly show it- 
self by increased non-circular motions. But this example is 
of course rather contrived, as for galaxies like NGC 2841, the 
center is usually so well-defined that uncertainties of 10" as 
modeled here are unlikely to occur within the THINGS sam- 
ple. 

We have shown that the results of a harmonic decomposi- 
tion of galaxies like NGC 2841 (i.e., having a steep rotation 
curve) are sensitive to offsets in the galaxies' center position. 
For NGC 2366, i.e., a dwarf galaxy showing close-to solid- 
body rotation, our results are mostly insensitive to small off- 
sets in the center position. The less well-defined center for 
NGC 2366 does therefore not affect our results of that galaxy. 

7. SUMMARY AND CONCLUSIONS 

We have analyzed H I velocity fields of 19 THINGS disk 
and dwarf galaxies at the currently best available spatial and 
spectral resolution. The observations, data reduction and anal- 
ysis of all these g alaxies were done i n an identical and homo- 
geneous manner dWalter et aT]|2008b Ide Blok et al.ll2008l) . We 
have determined the center positions of these 19 galaxies by 
comparing the kinematic center estimates with those from the 
radio continuum and/or NIR data. 

We have derived reliable center positions and show that 
most galaxies in our sample do not exhibit significant offsets 
between their kinematic and their optical centers. For 18 out 
of 19 galaxies from our sample, we have performed a har- 
monic decomposition of their velocity fields with the aim of 
quantifying the amount of non-circular motions and deriving 
a lower limit for the elongation of the potential. 

Our results show that in the (large) majority of our sam- 
ple galaxies, the effects of non-circular motions are small, 
irrespective of whether these non-circular motions are mea- 
sured in the inner few kpc or over the entire radial range. 
Within our sample, the average amplitude of the non-circular 
motions is A r = 4.8 ±4.0 kms" 1 for the inner 1 kpc and 
A r = 6.7 ±5.9 kms" 1 when averaged over the entire radial 
range. High non-circular motions in the inner parts of the 
galaxies are mostly found in barred, and/or luminous galaxies. 
The galaxies in our sample least dominated by baryonic mat- 
ter (i.e., the low luminosity galaxies) show the smallest non- 
circular motions in the sample, especially in the inner parts, 
thus not showing the streaming motions one would associate 
with elliptical distortions or tri-axial halos. 

The inner few kpc of the large majority of the galaxies in 
our sample contain non-circular motions of the order of ten 
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FIG. 17. — The results of the unconstrained (black dots) and constrain ed (gray crosses) harmonic decomposition of NGC 3198. For the constrained fit, the 
inclination and position angle were fixed to the values from de Blok et al. 12008). The layout of the figure is identical to that of Figs.|3]and|4] except that we do 
not show any weighted means in the panels on the right-hand side. The parameters that are plotted in this figure are defined in Sec. [3] The constant inclination in 
the inner pails of NGC 3198 cause a non-zero C3 term which results in slightly higher non-circular motions in the inner parts. 



percent of the local rotation velocity. Averaged over our sam- 
ple, the median amplitude of the non-circular motions in the 
inner 1 kpc is 8 ± 3 percent of the local rotation velocity. For 
three out of the four dwarf galaxies in our sample, the am- 
plitudes of the non-circular motions amount to (at most) ten 
percent of the local rotation velocity. 

Whether the non-circular motions are expressed in absolute 
terms or if they are normalized by the local rotation velocity 
does not affect the conclusions: the amount of non-circular 
motions detected in the majority of our galaxies is signifi- 
cantly smaller than what is expected from CDM simulations 
and they are f ar too small to "hide a cu sp in a core" as pro- 
posed, e.g., by Hav ashi & Navarrol (120061) . 

The average elongation of the gravitational potential and its 
scatter, both statistically corrected for the unknown viewing 
angle ipi, is (e pot ) = 0.017 ±0.020. This is significantly lower 
than the predictions from cosmological simulations. The indi- 
vidual elongation measurements show that the large majority 
of the galaxies in our sample have elongation measurements 
which are systematically below the CDM predictions. The 
elongation of all galaxies in our sample is consistent with a 
round potential, although some galaxies have large enough 
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FIG. 19. — Harmonic analysis of NGC 3198 using the hermite velocity field (black dots) and the intensity weighted mean velocity field (gray crosses). The 
layout of the figure is identical to that of Figs. [3] and \4\ except that we do not show any weighted means in the panels on the right-hand side. The differences 
between the derived quantities are small. The radial variation of the harmonic components seems, however, more pronounced in the hermite velocity field. 
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Fig. 20.— Like Fig.[T9lbut for DDO 154. 
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FIG. 21. — Harmonic analysis of NGC 2366 using our best center position (black) and a center position, offset by 10" along the major axis (gray). The layout 
of the figure is identical to that of Figs.fJ]and|4] except that we do not show any weighted means in the panels on the right-hand side. The difference between the 
two decompositions is marginal. The offset center causes the median amplitude of the non-circular motions (dashed line in the panel showing the distribution of 
A,-(r)) to increase fromA r ~ 3kms~' toA r ~ 3.5kms~', showing that a small offset along the major axis has no significant effect on the derived quantities in the 
case of NGC 2366. 
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FIG. 22. — Like Fig.ffT] but the harmonic analysis shown in gray is using a center position offset by 10" along the minor axis. Although being larger than in 
the case of the offset along the major axis, the difference between the two decompositions is still small. The offset center causes the median amplitude of the 
non-circular motions (dashed line in the panel showing the distribution of A r (r)) to increase from A, ~ 3kms~' to A r ~ 4 km s , showing that even a small offset 
along the minor axis does not influence the derived quantities in a significant way in the case of a dwarf galaxy like NGC 2366. 



20 



Trachternach et al. 






15 


cn 


10 


E 










5 


E 










2 



^,wt/-, r, 1 1 1 1 1 1 = 

200 400 600 
r (arcsec) 



30 
20 
10 




0.7 




200 400 600 
r (arcsec) 



FIG. 23. — Harmonic analysis of NGC 2841 using our best center position (black) and a center position, offset by 10" along the major axis (gray). The layout 
of the figure is identical to that of Figs. [3] and [4] except that we do not show any weighted means in the panels on the right-hand side. For NGC 2841, the 
amount of non-circular motions increases clearly when choosing an offset center. The largest differences are visible in the cj term. The median amplitude of 
the non-circular motions (dashed-line in the panel showing the distribution of A r (r)) increases from A,- ~ 7kms~' to A r ~ 12kms~' , showing that in the case of 
NGC 2841, an offset along the major axis will clearly show itself as an increase in the amplitudes of the harmonic components. 
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FIG. 24. — Like Fig. I231 but the harmonic analysis shown in gray is using a center position offset by 10" along the minor axis. For NGC 2841, the amount 
of non-circular motions increases clearly when choosing an offset center. The largest differences are visible in the $2 term. A r (dashed-line in the panel showing 
the distribution of A r (r)) increases from A r ~ 7kms~' to A,- ~ 25kms~', showing that in the case of NGC 2841, an offset along the minor axis will produce 
significantly different results than if the dynamical center position is chosen. 
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APPENDIX 

A. DESCRIPTION OF INDIVIDUAL GALAXIES 

In this section, we present for each galaxy in our sample a detailed description of (a) the derivation of the center position, and 
(b) the results of the harmonic decomposition. The galaxies are presented in order of increasing right ascension. The results are 
presente d in graphical for m in Figs. 12511421 Channel maps, optical images and moment maps for the galaxies in our sample are 
given in lWalter et al.l (120081) . The different center estimates are summarized in Table [2] where our adopted best center positions 
are shown in bold face. The results from the harmonic decomposition are summarized in Table [3] 

A.l. NGC 925 (Fig.\25) 

a ) Center estimates 

NGC 925 is a barred late-type spiral galaxy whose H I disk looks rather clumpy. The 3.6 /im IRAC image shows a weak central 
bar and two spiral arms in the outer region. NGC 925 lacks a clear center in the IRAC image as well as in the radio continuum. 
We used ELLFIT at varying isophote levels to derive a central position from the IRAC image. This center is offset from the 
pointing center by SX = -0"8 and SY = —2'.' 6 (i.e., towards the south-east) and was used as an initial estimate for a ROTCUR 
run with all parameters left free. Averaging the values for 5X, i5Y over the region not too heavily affected by spiral arms and 
asymmetries (r < 121", see top row of Fig. [25]) results in a center position offset from the pointing center by SX = 4." 4 ± 8."3 and 
SY = —l'.'S ±4"2. As this position is consistent with the center from the IRAC image (see row 2 of Fig. [25t. we adopt the former 
one as our best center for NGC 925. 

b) Harmonic expansion 

The fitted values for the inclination show a lot of scatter in the inner 100". This is partially caused by the near solid-body rotation 
in that region, which makes determining a kinematic inclination more difficult. Note that in the case of pure solid-body rotation, 
VROT and sin(/) are degenerate and cannot be fitted simultaneously. However, we seldom have to deal with pure solid-body 
rotation, but more usually with near solid-body rotation, which gives the algorithm some handle on the fitted values. Although 
the amplitudes of the individual components are not constant over the radial range, the median amplitudes of all harmonic orders 
are relatively small (A m < 5 kms" 1 when averaged over the entire galaxy). The median amplitudes are slightly higher if one 
averages over the inner 1 kpc only, which might be due to the stellar bar. A r is ~ 6.3 kms" 1 if averaged over the entire range and 
~ 9.5 kms" 1 if averaged over the inner 1 kpc only (see also Table [3). The elongation of the potential is unconstrained inwards 
of r w 100", mainly because of the large uncertainties for the derived inclination (which enters into the uncertainty in en 0t as a 
fourth power). Its weighted mean is consistent with a round potential ((e pot sin(2iy92)} = 0.000 ± 0.046, see again Table Bb. The 
median of the absolute residual velocity field is 3.0 kms" 1 , showing that a harmonic decomposition up to third order was able to 
capture most non-circular motions present in NGC 925. 

A.2. NGC 2366 ( Fig. \26$ 

a ) Center estimates 

NGC 2366 is classified as a dwarf irregular and belongs, like NGC 2403, to the M81 group. As for all the other dwarf galaxies 
in the sample presented here, the 3.6 /im IRAC image and the radio continuum map do not show a nuclear source. We derived the 
center of the emission in the IRAC image by using ELLFIT at varying isophote levels. The IRAC center is offset from the pointing 
center by 5X = 7"8 ± 0."6 and SY = -18"0 ± 1 "8. This center estimate was used for a ROTCUR run with all parameters left free. 
We averaged the derived values for <5X, SY for r < 200" (see row 1 in Fig. |26| |. but excluded the discrepant data points between 
70" < r < 100". The resulting dynamical center is offset from the pointing center by SX = 5."2 ± 3"5 and SY = —19 ."2 ± 7"7, 
and within the uncertainties consistent with the center position derived from the IRAC image (see row 2 of Fig. |26). Deriving a 
kinematic center by averaging over all data points out to r = 200" gives a consistent center estimate, though with a larger scatter 
(SX = 2"9 ± 6."5, SY = —18 ."5 ± 7"9). Therefore, we adopt the former kinematic estimate as our final center position. 

b ) Harmonic expansion 

NGC 2366 shows no clear signs of spiral structure. Its rotation curve is clearly dominated by near solid-body rotation, making 
it difficult to determine the rotation velocity and i simultaneously. At r ~ 100", the systemic velocity co rises from ~ 100 kms" 1 
to ~ 105 kms" 1 . At the same radius, the amplitude of 52 "jumps" from ~ kms" 1 to ~ 4 kms" 1 and decreases then to -4 kms" 1 
in the outskirts of NGC 2366. The distribution of A r (r) shows that the largest amount of non-circular motions is to be found in 
the outer parts of NGC 2366, and that the inner 1 kpc exhibit only minor non-circular motions. The weighted mean elongation is 
(e P ot sin(2< / 92)} = 0.003 ±0.067 and the median of the absolute residual velocity field is 2.4 kms" 1 , once again showing that the 
galaxy has a round potential and that there are no large non-circular motions present in the m > 3 terms. 

A.3. NGC 2403 ( Fig. [27} 

a ) Center estimates 

NGC 2403 belongs to the M81 group, and is a late-type Sc spiral. Its 3.6 /im IRAC image shows multiple spiral arms and a 
bright central component. However, neither the IRAC image, nor the radio continuum map show a nuclear source. Therefore, 
we used 3.6 /im isophote fits to determine center coordinates and used the resulting center position (<5X=-4"7, SY=— 13"6) as 
an initial estimate for an unconstrained ROTCUR run. As can be seen in the first row of Fig. [27] SX and SY stay fairly constant 
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over a large part of the galaxy. We determine the dynamical center by averaging the values for r < 220", thus excluding the 
radii where SX changes considerably, probably due to the spiral arm located at that radius. The offset from the pointing center 
is SX = -5"2 ±5."2 to the east and SY = -12"8 ±4."1 to the south. Averaging SX and <5Y over a larger part of the galaxy gives 
comparable results but with a larger scatter. Given the agreement between the dynamical center estimate and the center as derived 
using the IRAC image, we adopt the r < 220" kinematic center estimate as our best center position. 

b ) Harmonic expansion 

The radial variation of the inclination is only ~ 5°, except at a few inner radii. The PA shows a steep inner rise, a fiat outer 
part and a dip at r « 200", which coincides with a high amplitude in the C2 and S2 components. Inspection of the integrated 
H I map shows a spiral arm crossing at this radius. The s\ and S3 components also show the characteristic wiggles caused by 
spiral arms. Additionally, S3 has a small offset from zero, which through Eq. [3]can indicate a slightly elongated potential. The 
average elongation of NGC 2403, however, is (e pot sin^^)} = -0.022 ±0.025, and therefore consistent with a round potential. 
The distribution of A m shows that the individual harmonic orders contribute only ~ 2 percent to V tot , indicating that non-circular 
motions play a minor role in NGC 2403. The latter can also be seen in the distribution of A r (r). 

A.4. NGC 2841 (Fig.\28b 

a ) Center estimates 

NGC 2841 is an early-type spiral galaxy. Its H I distribution has a central hole and the IRAC 3.6 fim image shows a distinctive 
bulge and a flocculent spiral structure. The radio continuum map of NGC 2841 shows a strong nuclear source coinciding within 
0"2 with the position of the central source in the 3.6 /im IRAC image. We used the position of the radio continuum source 
(<5X=-6."3, (5Y=0."5) as input for ROTCUR and made a fit with all parameters left free. We averaged the values for 5X and SY 
between 75" < r < 220" (see first row in Fig.l28l. The lower limit excludes the two innermost rings, as these are sparsely filled. 
The upper limit restricts our averaging to a region which is unaffected by spiral arms or by the warp. The dynamical center 
derived in such a way (SX = —6 ."5 ± 0"6, SY = 0"4 ± 0."9) deviates by less than 0."2 from the other center estimates. Therefore, 
we choose the radio continuum center as the best center position. 

b) Harmonic expansion 

The inner 200"-250" of NGC 2841 are in unperturbed circular rotation and show only small harmonic terms. For the inner 1 
kpc, we have no data as the H I surface density in this region falls below the 3cr detection limit imposed on the velocity fields. 
Most of the parameters show a large change at r ks 350" -400": the systemic velocity drops by about 20 kms" 1 , and most of the 
harmonic components show their highest amplitude at these radii. Looking at the total H I intensity map, one can see that this 
radius coincides with the location of a strong spiral arm. The dip in the fitted inclination at those radii causes a non-zero C3 term, 
induced by the spiral arm. The radial variation of the PA and inclination indicates that the outer disk of NGC 2841 is warped. All 
harmonic components are similar in amplitude, having median amplitudes of 3 < A„, < 4 km s" 1 , or < 2 percent of V tot . A r (r) is 
small in the inner parts, but > lOkms" 1 in the outer parts ofNGC 2841. Its median value is A r ~ 7 kms" 1 . The elongation of the 
potential is fairly constant over the radius and its weighted average is consistent with a round potential (see Tabled. 

A.5. NGC 2903 (Fig.\2§ 

a ) Center estimates 

NGC 2903 is a barred galaxy with tightly wound spiral arms. Its dominant bar can be seen both in the IRAC image and in the 
radio continuum map. The radio continuum additionally shows a strong nuclear source, coinciding within 1" with the center of 
the bulge in the 3.6 /im image. The central position as given by the radio continuum (SX=-1'.'2, SY=2'.'4 ) was used to make a 
ROTCUR run with all parameters left free. To derive the kinematic center, we averaged the SX and SY values over the stable parts 
of the fit, i.e., for 100" < r < 300" (see top row in Fig. |29). The resulting kinematic center (SX = 0"4 ± 2."4, SY = 0"6 ± 2"2) 
agrees within the uncertainties and within one beam with both the IRAC and the radio continuum center (cf. second row in 
Fig. 129). We therefore adopt the center position as derived from the radio continuum as our final central position. 

b ) Harmonic expansion 

NGC 2903 shows high amplitudes in all harmonic components in the inner 125". These non-circular motions are presumably 
caused by the bar in NGC 2903. The innermost 4-5 data points show the largest deviations from the general trends, most clearly 
visible in the distributions of the PA, inclination and C(>. 

The pronounced shape in the distribution of so could be caused by an m = 1 or m = 3 term in the gravitational potential. Given 
that NGC 2903 is not kinematically lopsided (there is, e.g., no significant difference in the rotation curves of the approaching 
and receding sides ), the contribution from an m = 1 term is likely to be small. Inspection of the total H I maps presented in 
Walte ret al.1 d200 8) shows that a spiral arm is located at the same radii as the pronounced variation in S2. The wiggles in the radial 
distributions of si and S3 at r ~ 100" and r ~ 350" also coincide with the locations of spiral arms. The median amplitudes of 
the individual harmonic components are small when averaged over the entire galaxy (A m < 4 kms" 1 , or 2-3 percent of V tot ), but 
increase towards the center of the galaxy. Within the inner kpc, the m = 2 component of A m has a median value of about 1 4 km s" 1 . 
These high non-circular motions, which can also be seen in the distribution of A r (r), are likely to be caused by the strong bar, 
and the associated streaming motions. Like some of the c,, s,, and thus also A r (r), the elongation of the potential shows some 
variations with radius. Its mean value is nevertheless consistent with a round potential (see Table|3]l. 
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A.6. NGC 2976 (Fig.\M 

a ) Center estimates 

The 3.6 /im image of NGC 2976 shows no sign of a bar or spiral arms. Nevertheless, the image shows two enhanced star 
forming regions at either end of the disk which coincide with density enhancements in the total H I map. Additionally, the IRAC 
image contains a nuclear source in the central parts of the galaxy, presumably a nuclear star cluster, which has, however, no 
counterpart in the radio continuum. To test whether the position of this central source agrees with the dynamical center, we 
perform an unconstrained ROTCUR run and use the position of the nuclear star cluster (<5X= 0."2, SY= 0."1) as an initial center 
estimate. By averaging SX, SY over r < 50" (see Fig. [30] top row), we derive a dynamical center which is offset from the 
pointing center by SX = 2."2 ± 3"0 and SY = 0."8 ± 1 "4 towards the north-west. Averaging over the entire galaxy gives consistent 
results (SX = 2."5 ± 3."0, SY = l."6± 3"9). As the dynamical center and the position of the nuclear star cluster agree within the 
uncertainties and to within one beam (cf. second row in Fig. [30b . we adopt the position of the nuclear star cluster as the best 
center position of NGC 2976. 

b) Harmonic expansion 

The harmonic components of the velocity field of NGC 2976 show a regular behavior. Both s\ and 53 scatter within ~ 1 kms" 1 
around zero, indicating that the elongation of the gravitational potential is very small. Its weighted mean is consistent with zero 
((e pot sin(2^ 2 )} = -0.010 ± 0.018). 

The C2 component is slightly offset from zero and shows sine-like variat ions with radius, especially for r > 80", where the two 
H I density enhancements are located. According to lSchoenmakersId 19991) . a non-zero m = 2 harmonic component in the velocity 
field can be caused by an m = 1 or m = 3 term in the potential. NGC 2976 does not se em to be kinematic ally lopsided (cf. the 
negligible differences in the rotation curves of the approaching and receding sides in lde Blok et al.l |2 008i) . Therefore, there is 
most likely only a small contribution of an m = 1 term in the potential to the m = 2 harmonic term in the velocity field, and the 
significant contribution should come from an m = 3 term in the potential. 

Beyond r ss 80", the PA stays constant. The median amplitudes of all harmonic component are small (A m < 3 percent V tot ), 
irrespective of whether one averages over the entire galaxy or only over the inner 1 kpc. The distribution of A r (r) is also small 
for most radii. Only a few data points at extreme radii show large amplitudes, but these are associated with large uncertainties. 
The mean elongation of the potential is, despite being unconstrained in the inner 30", again consistent with a round potential (see 
Table©. 

A.7. NGC 3031 (Fig.\JB 

a ) Center estimates 

NGC 3031 (better known as M 81) is a grand design spiral. Its two well-defined spiral arms are easily visible in the IRAC 
3.6 [im image as well as in the H I intensity map. NGC 3031 contains a strong central radio source. The 3.6 /im image also shows 
a well-defined central component with a central minimum which coincides with the position of the central radio-continuum 
source (SX= -202"7, SY= -371. "9). We used the position of the radio continuum source as input for a ROTCUR run with all 
parameters left free. As can be seen in the top row of Fig. [31] SX and <5 Y show a clear break at r ~ 400", probably because of the 
prominent spiral arm which is located at approximately that radius. We therefore estimated the dynamical center by averaging 
SX and SY over 200" < r < 385", thus omitting the innermost data points whose tilted-rings are only sparsely filled and have 
correspondingly large uncertainties. The resulting dynamical center [SX = -205 "1 ± 3"1, SY = -374 "6 ± 3. "9) coincides within 
the uncertainties and within less than one beam with the other center estimates (second row of Fig.[3TI) and we therefore adopt 
the position of the radio continuum source as the best center position of M 81. 

b ) Harmonic expansion 

The PA, co, and the inclination of M 8 1 rise steadily beyond r ss 400", indicating that the outer disk is warped or disturbed, 
which is not unexpected given that M 81 interacts with M 82 and NGC 3077. The radii with the highest gradient in s\ and S3 
(r Ri 350" -400") correspond with the location of a prominent spiral arm. Between 400" < r < 700", and $2 change from 
values around +5 kms" 1 to < -10 kms" 1 . Note that if we were to perform a harmonic decomposition with an unconstrained 
center position, the center would vary at those radii in order to minimize the ci and S2 term, thus underestimating the amount of 
non-circular motions present in M 8 1 . However, as we keep the center position fixed, we are able to quantify and detect these 
non-circular motions. The median amplitudes A„, of the m = 2 and m = 3 component are roughly equally high (~ 5 kms" 1 ), 
whereas that of the m = 1 component is slightly smaller (~ 2 kms" 1 ). As M81 has a central H I minimum, we have no data for 
the inner 1 kpc. The radial distribution of the non-circular motions varies between 3 < A r (r) < 15 kms" 1 and has a median of 
A, ~ 9 kms" 1 (<5 percent of V tot ), making the amount of non-circular motions quite small, despite the prominent spiral arms. The 
elongation of the potential shows distinct radial variation which is due to the spiral arms. Its weighted average is again consistent 
with zero (see Table©. 

A.8. NGC 3198 (Figs.UM 

a) Center estimates 

NGC 3 198 is classified as a SBc spiral. Its 3.6 /im IRAC image shows two well-defined spiral arms, emanating from a prominent 
bulge. The central component has a nuclear point source embedded, which has a counterpart in the radio continuum. The center 
estimates from the IRAC image and the continuum map agree to within 1". We make a ROTCUR run with the position of the radio 
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continuum center (<5X= —I'M, SY= -0"1) as an initial estimate. Fig.Q]shows that the outer parts (r > 150") of NGC 3198 are 
strongly affected by the spiral arms. For the derivation of the dynamical center, we therefore restrict the averaging of <5X and SY 
to radii with r < 100". The derived dynamical center (SX = — 1."4± l'.'l, SY = —0'.'2 ±2."1) agrees well with the other estimates 
(see Fig. |2). We therefore adopt the position of the continuum source as the best center position for NGC 3 198. 

b) Harmonic expansion 

The PA of NGC 3198 shows a steep increase within the inner 200", and then declines slowly. The inclination varies in the 
inner parts in a range of about five degrees, but shows a steady increase beyond r ~ 450", indicating that the outer disk is 
warped. The C3 term is small, meaning that the fitted inclination is close to the intrinsic inclination of the disk. Although there 
is no global offset from zero for c 2 and s 2 , they show small deviations from z ero for some radii w hich are coinciding with the 
locations of spiral arms in the 3.6 fim image or the total H I map presented in lWalter et alj d2008l) . The distributions of s\ and 
S3 are best described as wiggles caused by spiral arms on top of a small offset. The global elongation of the potential is with 
(e po t sin(2y> 2 )) = 0.017 ± 0.020 consistent with zero. 

The median amplitudes of the individual harmonic components are similar in amplitude, and of order A,„ ~ 2-3 kms" 1 , or < 
2 percent of V tot . The distribution of A r (r) shows that the amplitude of the non-circular motions is < 8 kms" 1 for most radii. Its 
median is A r ~ 4.5 kms" 1 when averaged over the entire radial range, and A r ~ 1 .5 kms" 1 for the inner 1 kpc. The median value 
of the absolute residual velocity field is ~ 2.6 kms" 1 , showing that a harmonic decomposition up to third order was capable of 
capturing most non-circular motions. 

A.9. IC2574(Fig.\32b 

a) Center estimates 

IC 2574 has neither a central radio continuum source, nor a clear nuclear source in the 3.6 IRAC image. The dynamics of this 
galaxy show clear evidence of random non-circular motions, showing themselves as kinks in the iso-velocity contours . For the 
cente r determination, we make use of a velocity field where most non-circular features were removed (as p resented in lOh et ail 
2008). We make a ROTCUR run with all parameters free, using the center position derived by lOh et ail d2008l) as an initial center 
estimate. The results of the ROTCUR fit are shown in row 1 of Fig. [32] It can be seen that the values for <5X and SY do not 
change much over a large part of the galaxy. Averaging SX, SY over r < 400" results in a dynamical center which is offset by 
SX = 8"7 ± 14."9 and SY = 18. "6 ± 10"4 towards the north-west of the pointing c enter. Thi s corre sponds to a 2 ooo = 10''28 m 27.5 s , 
(52000 = +68 °24'58"7 and is in excellent agreement with the estimate derived bv lOhet all J2008) (a 2 ooo = 10 /! 28 m 27.7' 5 , 62000 = 
+68°24'59."4). 

b) Harmonic expansion 

Although we made use of the bulk velocity field from lOh et al.1 (120081) for the center estimate, we use the hermite velocity field 
for the harmonic decomposition, as we aim to quantify non-circular motions. The results of the harmonic decomposition of the 
hermite velocity field of IC 2574 are shown in the third row of Fig. [32] 

At a radius of r ~ 250", the PA and inclination drop significantly, and the harmonic components jump to more extreme 
values. The distribution of A r {r) has its maximum at this radius, and c\ also shows a clear break. The fe ature which is mos t 
likely responsible for the major non-circular motions is the supergiant shell in the north-east of IC 2574 (see lWalter et al.l fl998). 
Looking only at the inner 1 kpc, the amplitudes of the individual harmonic components are all < 1 kms" 1 . If one measures 
the amplitudes over the entire radial range of the galaxy, the m = 2 component is the dominant one, mainly because of the large 
amplitudes in c 2 and s 2 at r ~ 250". 

The median amplitude of the quadratically added non-circular motions is A r ~ 3.8 kms" 1 when averaged over the entire radial 
range, and A r ~ 1.4 kms" 1 for the in ner 1 kpc. This shows that although IC 2574 contain a significant amount of chaotic non- 
circular motions (see lOh et al.ll2.Q08b . the effect of systematic (potential induced) non-circular motions is small. For some parts 
of the galaxy, the elongation of the potential is unconstrained. Its weighted mean is consistent with a round potential, although it 
has a large uncertainty associated ((e pot sin(2i^ 2 )} = 0.015 ±0.045). 

A. 10. NGC 3521 (Fig.\33) 

a ) Center estimates 

NGC 3521 appears in the 3.6 /im IRAC image as well as in the H I map as a multi-armed, flocculent spiral galaxy. The IRAC 
image shows an inner disk, in which a strong nuclear point source (without a counterpart in the radio continuum) is embedded. 
We use the position of the nuclear point source (SX= -8 "8, SY= 5 "7) as input for a ROTCUR fit with all parameters left free. 
The top row of Fig. [33] shows that the values of SX and SY stay constant for radii smaller th an 220". Beyond that radius, the 
rotation curves for the approaching and receding sides start to differ (see Ide Blok et al.ll2008l) and this most likely also affects 
the derived central positions. In order to derive a dynamical center, we average SX and SY for r < 220". The resulting center 
(SX = -8"9 ± 0."8, SY = 6."6 ± 1 "4) agrees with the central position found in the IRAC image (see second row of Fig. l33l and 
we adopt the IRAC position as the best center estimate. 

b) Harmonic expansion 

As can be seen, e.g., in the radial distribution of PA, i, and A r (r), the inner few data points clearly deviate from the rest. Note 
that the C3 component is higher for these points, meaning that their inclination could not be fitted correctly in an unconstrained fit. 
Some of these inner data points show large non-circular motions, while others (e.g., the inner two points) do not. This becomes 
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even clearer if one looks at the distribution of A m . Although the median of the m = 3 component in the inner 1 kpc is small, 
its upper quartile is rather large, as the inner 1 kpc contain only three data points, one of which has a large amplitude. The Si 
and S3 components are both slightly offset from zero and superimposed on a wiggle at r « 200". The wiggle coincides with a 
density enhancement in the total H I intensity map (in this case a spiral arm). However, as si is offset to negative values, and S3 
to positive ones, their offsets almost cancel out for the derivation of the elongation of the potential (cf. Eq.0, whose weighted 
mean is within the uncertainties consistent with zero ((e po t sin(2ip2)} = 0.017 ± 0.019). 

The effect of the spiral arms is clearly visible in the C2 and S2 components. NGC 3521 has the highest median value in the 
absolute residual velocity field of our sample (4.5 kms ). 

A. 11. NGC 3621 (Fig.\34b 

a ) Center estimates 

The 3.6 /im IRAC image of NGC 3621 is dominated by a flocculent spiral structure, which contains a central point source. 
The prominent star forming region in the south-west is connected to the center of the galaxy by a spiral arm which is clearly 
visible in the total H I map and the IRAC image. While this star forming region has a counterpart in the radio continuum, the 
IRAC nuclear point source has not. We use the coordinates of the nuclear point source in the IRAC image (SX= -4 ."2, 8Y= 
-1"1) as an initial estimate for a ROTCUR fit with all parameters left free. The results (top row of Fig.l34l show that the derived 
central positions are fairly constant for radii smaller than 450" and we therefore derive the dynamical center by averaging over all 
points with r < 450", with the additional exclusion of the innermost data point. The resulting dynamical center (5X = -6"0±2"0, 
SY = 1"4±6."5) coincides with the coordinates of the point source both within the uncertainties and the size of the beam (Fig. [34] 
second row). Therefore, we adopt the position of the central point source as our best center estimate for NGC 362 1 . 

b ) Harmonic expansion 

The data points beyond r = 600" are not used for our analysis due to their sparsely filled tilted-rings. The PA of NGC 3621 is, 
except for the innermost two data points, constant to within a few degrees. Apart from a few data points at extreme radii, the 
amplitude of s\ is almost negligible. C3 is close to zero, except for r « 300". This position also coincides with the radius of the 
highest amplitude in C2, and si- Inspection of the total H I intensity map shows that there is a ring-like H I density enhancement 
at which a spiral arm emerges in the south-east. The S3 component shows clear wiggles, indicating spiral arms, and a steep rise 
beyond r w 500", where the tilted-rings are less filled. The elongation of the potential again shows the kinematic signature of 
the spiral arms. Its weighted average is consistent with zero (see TableQ. The median amplitudes for each individual harmonic 
component are fairly small (A,„ < 3 km s" 1 , or < 2 percent of V tot ). The distribution of A r (r) is < 5 km s" 1 for the majority of radii 
and has a median value of A r ^3.4 kms" 1 . 

A. 12. NGC 3627 (Fig.\35$ 

a ) Center estimates 

The 3.6 /im IRAC image shows that NGC 3627 is a barred galaxy with two asymmetric spiral arms, which are also clearly 
visible in the total H I map. The central parts of NGC 3627 contain a nuclear point source which is visible in the IRAC image 
as well as in the radio continuum. Using the position of the point source {SX = 0"1, SY = 0"1) as input for an unconstrained 
ROTCUR fit results in a slightly different center estimate (see row 1 and 2 of Fig. 135V We exclude some data points at extreme 
radii because of their large uncertainties, and average <5X, SY for 45" < r < 160", resulting in a dynamical center which is offset 
from the pointing center by 5X = -4."5 ± 3"1 and SY = -6"8 ±4."8. These values deviate from the radio continuum source by 
approximately 1.5(7. Given the asymmetric appearance of NGC 3627, it is to be expected that the results from ROTCUR are 
somehow affected by it. However, the difference between the kinematic center and the center from the radio continuum is only 
one beam. We therefore use the coordinates of the central radio continuum source as the best center position for NGC 3627. 

b) Harmonic expansion 

NGC 3627 has a central hole in the H I distribution. Therefore, we see only the flat part of the rotation curve and have no 
data within the inner kpc. The systemic velocity changes from ^730 kms -1 at r < 60" to ^705 kms -1 for r > 90". All 
harmonic compone nts (C2, C3, Si, S2, S3) show high amplitudes. The galaxy is kinematically and morphologically lopsided (cf., 
Ide Blok eT al. 2008), which is also indicated by the non-zero values for cj and S2. Both the S\ and the S3 component show a "dip" 
at r 70", followed by a steep rise beyond that radius. In the case of S3, the rise is followed by a second decline. The radius of 
that first "dip" coincides with the location of the large arm which extends to the south of NGC 3627. The median amplitudes of 
the individual harmonic components are the highest in our sample, reaching from A m ~ 10 kms" 1 for m = 3 to A m ~ 17 kms" 1 
for m = 2. The distribution of A r (r) shows that the amplitude of non-circular motions is large at all radii. The median amplitude 
is A r ~ 28.5 kms" 1 , or up to 14.7 percent of V tot . Given the morphology of NGC 3627, this is not unexpected. Despite the 
high non-circular motions present here, the weighted mean elongation of the potential is still small, and within its uncertainty 
consistent with a round potential (cf. Table [3), although there are large errors associated with the data points in the outer parts 
and in the region at r ss 70". 

A.13. NGC 4736 (Fig.\M 

a ) Center estimates 

The 3.6 ^m IRAC image of NGC 4736 is dominated by a ring of star formation which is also prominently visible in the radio 
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continuum and the total H I map. The center is well-defined both in the 3.6 /im image and in the radio continuum. The former 
shows a small central minimum in the flux which coincides within l."3 with the central radio-continuum maximum. We use the 
position of the radio continuum source [SX = 0"2, SY = — 0."4) for a ROTCUR run with all parameters left free. The fitted center 
positions are stable within the inner 100", but have an increasingly large scatter towards the outer parts, most likely due to the 
marginally filled tilted-rings beyond r ~ 100", and also because of the low inclination of NGC 4736. We therefore average the 
estimates for <5X and <5Y over the inner 100", resulting in a dynamical center (SX = 0"5 ± 1 "8, SY = 0."6 ± 2."0) whose deviation 
from the other center estimates slightly exceeds the uncertainties (1.2er deviation for <5Y). Though, to put things into perspective, 
this is far smaller than one beam (see second row of Fig. [36). Thus, we will use the position of the radio continuum source as the 
best center estimate for NGC 4736. 

b ) Harmonic expansion 

The star forming ring at r ~ 80" causes a lot of confusion in the inner parts of NGC 4736. Here the unconstrained fit by RES WRI 
produces an apparent rotation velocity which rises to twice the value in the outer parts. The other harmonic components also 
show high amplitudes. The median amplitudes of each harmonic component (averaged over the inner 1 kpc) are negligible for 
the m = 1 and m = 3 component (A m ~ 1 kms" 1 ), but much higher for the m = 2 component (A m ~ 7 kms -1 ). The amplitude of 
the quadratically added non-circular motions, A r (r), is ~ 10 km s" 1 for most radii. 

NGC 4736 is both kinematically and morphologically lo psided, which can b e seen not only in the differences between the 
rotation curves of the approaching and receding sides (see Ide Blok et al.ll2008l) . but also in the high amplitudes in the C2 and 
52 terms of our harmonic expansion. The highest amplitudes in the harmonic terms are found for radii between r w 200" and 
r ss 300", i.e., those radii where the large northern spiral arm is located. The weighted average of the elongation of the potential 
is (e pot sin(2p 2 )) = -0.055 ±0.149. 

A. 14. DDO 154 (Fig. [23* 

a) Center estimates 

DDO 154 is a quiescent, gas-rich dwarf galaxy with a warp in the outer parts. It does not have a compact, well-defined center 
in the 3.6 /im IRAC image, or in the radio continuum. We attempted to derive a center position by fitting ellipses at a few 
representative intensity levels to the IRAC image, and the total H I map. The resulting central positions showed a variation of up 
to 15", demonstrating that deriving an unambiguous photometric center for DDO 154 is not straightforward. We used the center 
which was derived using the total H I map as a initial estimate for a ROTCUR fit with all parameters left free. The central positions 
as derived by ROTCUR are fairly stable over the radial range of the galaxy and the kinematic center was derived by averaging 
the SX, SY values over 100" < r < 260", thus omitting the slightly more unstable inner and outer regions. The resulting center 
position is given in Table [2] Including the inner data points gives a similar center position, but with two times larger scatter. We 
therefore adopt the dynamical center as averaged over 100" < r < 260" as the best center position of DDO 154. 

b ) Harmonic expansion 

The C3 term is close to zero for the entire radial range, normally indicating that the inclination could be well-determined. 
However, due to the nearly solid-body rotation in the innermost parts of DDO 154, the inclination is not very well-constrained 
in this regio n, as can be seen in the large scatter of the individual data points for r < 70". This is studied in more detail in 
Section 16.41 The component shows a maximum at around r w 80" which coincides with a minimum in the 52 component 
and a change in the systemic velocity cq, providing evidence for a small kinematic lopsidedness as supported by the variation 
of the center positi on at these radii and the differences between the approaching and receding sides of the rotation curve (see 
Ide Blok eT al. 2008). However, the effects described above are small. Despite the small-scale structure present in the velocity 
field, the amplitudes of all harmonic components — both for the entire galaxy and for the innermost 1 kpc — are < 1 kms -1 , 
or ~ 2 percent of V tot . The distribution of A r (r) shows that the amplitudes of non-circular motions in DDO 154 are small for all 
radii. The mean elongation of the potential is small, (e pot sin(2iy92)) = 0.024 ± 0.033, and consistent with zero. The median of the 
absolute residual velocity field is 1.2 kms" 1 (the smallest in our sample). 

A.15. NGC 4826 (Fig.\38} 

a) Center estimates 

NGC 4826 shows no sign of spiral structure in the 3.6 /im IRAC image, but has two counter-rotating gas disks (Bra un et al.1 
1994). The center of NGC 4826 is well-defined by the central, compact source visible in the radio continuum, and the IRAC 
3.6 pm image, whose coordinates agree to within 1". We use the position of the radio continuum source (SX = 2."4, SY = 8."5) as 
an initial estimate for an unconstrained ROTCUR fit. We derive the dynamical center by averaging the <5X, SY values inwards of 
60", i.e., over the radial range of the inner disk. The kinematic center derived in such a way (5X = 2"2 ± 1 '.'6, SY = 7."5 ± 0."8) 
agrees well with the other center estimates and we therefore assume the position of the radio continuum source to be the best 
center for NGC 4826. 

b ) Harmonic expansion 

Because of the sparsely filled tilted-rings, and the additional complication of the two counter-rotating disks, we have not at- 
tempted to derive a harmonic decomposition for this galaxy. 
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A. 16. NGC 5055 (Fig. \3<fy 

a ) Center estimates 

NGC 5055 is classified as an Sbc galaxy and has a flocculent structure in the 3.6 /im IRAC image. The IRAC image shows a 
compact nucleus with a well-defined center, whose position coincides to within 1" with a faint central radio continuum source. 
We use the position of the latter (SX = 8f 8, SY = 0."4) as the input center for a ROTCUR run with all parameter left free. 

We derive the dynamical center by averaging the SX, SY values for radii smaller than 450" to exclude the parts with sparsely 
filled tilted-rings beyond r ~ 450". This center (SX = 8"2±4."8, SY = 0"4±2"0) coincides within the uncertainties and to within 
a beam size with the other center estimates. We therefore adopt the position of the radio continuum source as the best estimate 
for the central position of NGC 5055. 

The H I distribution and kinematics of NGC 5055 were recently analyzed by Battaglia et al.1 (120061) . They define a dynamical 
center coinciding with their optical center at 0:2000 = 13*15 m 49.25*, £2000 = +42°01'49."3, which is 4" to the north of our choice. 
As our choice was also guided by the position of the central source in the 3.6 /im image, this would imp ly that the respective 
photometric centers must be shift ed by 4" wi t h resp ect to each other. The optical center position listed in Batta glia et al.l (120061) 
is ultimately based on a listing in Maoz et ail (1 1996b of the position of the central UV source in NGC 5055 determined using an 
HST FOC image. On examination of the relevant image in the Hubble Space Telescope archive, we find the central source to be 
located at 0:2000 = 13 /! 15'"49.3'\ <52ooo = +42°01'46"2, which is in much closer agreement (Of 8) with our position. 

b) Harmonic expansion 

For the harmonic analysis of NGC 5055, a distinction between the inner and the outer parts has to be made. NGC 5055 has 
a well-defined H I disk which extends to about 450". Beyond that radius, the H I column density is low er and generally fall s 
below the 3 a column density limit which we imposed during the construction of the velocity fields (see de Blok et al. 2008). 
Therefore, the uncertainties increase rapidly and we considered only data points with radii r < 450" for the radial averaging 
of the parameters. The median amplitudes of the individual harmonic components (averaged over r < 450") are all fairly low 
(A m < 3 kms" 1 or ~ 2 percent of V tot ). The median amplitudes averaged over the inner 1 kpc contain only three data points and 
are also low for the m - 1 and m = 3 components, but quite high (~ 8 kms" 1 ) for the m = 2 component. The latter has a large 
scatter associated with it, mainly because of the innermost data point which shows high non-circular motions, but is derived from 
an only partially filled tilted-ring. The global elongation of the potential is well-constrained in the inner part of NGC 5055. The 
weighted mean elongation within 450" is again consistent with a round potential (see TableO. 

A. 17. NGC 6946 (Fig. goj 

a ) Center estimates 

NGC 6946 is a late-type spiral galaxy with well-defined multiple spiral arms. The IRAC 3.6 /im image shows a point source 
in the center of the galaxy which coincides to better than 1" with the central radio continuum source. We use the position of 
the latter (SX = -If 8, SY = 14."5) as input for a ROTCUR fit with all parameters left free and determine the dynamical center by 
averaging the SX, SY values for 22" < r < 290", thus excluding the two innermost points which are affected by sparsely filled 
tilted-rings. The resulting dynamical center (SX = -3 "3 ±4."4, SY = 1 1 "8 ± 5 "8) coincides within the uncertainties and to within 
one beam with the other center estimates (see Fig. |40] row 2). The large uncertainties in the kinematic estimate are caused by 
the relatively low inclination of NGC 6946, which makes an unconstrained tilted-ring fit of the velocity field more difficult. We 
adopt the position of the central radio continuum source as the best center for NGC 6946. 

b ) Harmonic expansion 

NGC 6946 is the galaxy in our sample with the lowest inclination. This introduces larger uncertainties in the rotation curve 
and the harmonic decomposition. For our analysis of the median amplitudes, and for the derivation of the average elongation, 
we regard only the inner 420", i.e., the region where the unconstrained apparent circular velocity c\ does not rise to a value two 
times that in the flat part. This radius also coincides with the locatio n from which, whe n moving outwards, the approaching and 
receding sides of the velocity field start to differ substantially (see lde Blok et alj |2008). The median amplitudes averaged over 
r < 420" are small for each harmonic component (A m < 5 km s" 1 ). Due to the deficiency of H I in the center of NGC 6946, we 
can not study the kinematics within the inner 1 kpc. Again averaging only over r < 420", the median of the quadratically added 
amplitudes of all non-circular components is A r w 7 kms -1 . The steep rise in 52 and 53 at r w 120" is probably caused by the 
spiral arm which crosses at that radius. Due to the low inclination of NGC 6946, and the associated large error in the derivation 
of the fitted parameters, the elongation of the potential is not well-constrained for r < 150" and r > 420". Its weighted mean 
(again averaged over r < 420"), however, is consistent with a round potential. 

A. 18. NGC 7331 (Fig.EJ} 

a ) Center estimates 

The 3.6 /im IRAC image of NGC 7331 shows prominent spiral arms, and a well-defined nucleus which has no counterpart in the 
radio continuum. We use the center from the 3.6 /im image (SX = Of 1, SY = Of 3) as input for an unconstrained fit with ROTCUR. 
Ignoring the two innermost points as their tilted-rings are only sparsely filled, we average the values for SX and SY over the 
stable region with r < 185". This results in a dynamical center of SX = Of ± 0f6, SY = -1 f 8 ±2."2, which agrees within the 
uncertainties and to within one beam with the IRAC center. Therefore, we adopt the IRAC center as our best center position for 
NGC 7331. 
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b) Harmonic expansion 

Although PA and inclination do not change much over the radial range (~4-6 degrees), they both increase outwards, indicating 
a warping of the disk. The small cj values denote that the inclination could be well-determined for most radii. The s\ and 
53 values both show a wiggle at r « 150", coinciding with the location of a spiral arm visible in the total intensity H I map 
and in the IRAC 3.6 /im image. The variation of cq, ci, and 52 shows that NGC 7331 is kinematica lly lopsided, as sho wn by 
the differences in the rotation curves of the receding and approaching sides of the velocity field (cf.. Ide Blok et alj|2008l) . The 
median amplitudes of the individual harmonic components are all small (A m < 5 kms" 1 ). Because of the central H I deficiency 
of NGC 7331, there is no kinematic data for the inner 1 kpc. The distribution of A,(r) varies between A r (r) ~ 2 kms" 1 and 
A r {r) ~ 10 kms" 1 and has a median of A r ~ 6 kms" 1 . The elongation of the potential is well-determined and close to zero for 
most radii. The wiggle at r « 150" is again an indication of the aforementioned spiral arm. The weighted mean elongation is 
{e po t sin(2<£>2)) = -0.003 ± 0.017, and therefore consistent with a round potential. 

A. 19. NGC 7793 (Fig.]42$ 

a) Center estimates 

NGC 7793 is a fiocculent spiral whose 3.6 /im IRAC image shows a well-defined central source without a counterpart in the radio 
continuum. We adopt the position of this source (8X = — 9."1, 5Y = -0"9) as the initial estimate for an unconstrained ROTCUR fit, 
and average the <5X, <5Y values over the stable part of the galaxy (r < 180"). The resulting kinematic center (SX = -9."9 ±2."6, 
5Y = l."8 ±2"1) agrees reasonably well with the IRAC center (1.5er deviation). However, the deviation between these two center 
estimates is ~ 6 times smaller than the size of one beam and we therefore use the IRAC center as our best center position. 

b) Harmonic expansion 

The PA and the inclination vary continuously with radius, indicating that the disk of NGC 7793 might be warped. The systemic 
velocity, co, stays fairly constant over radius. The offset from zero in the $i component might indicate a slight kin ematic lopsid- 
edness , which is also visible in the difference between the receding and approaching sides of the rotation curve in de Blok et al. 
(2008). The 53 term is significantly larger than the S\ term, which is typical for galaxies with low inclination dSchoenmaken 
1999). In fact, the 53 component seems to be offset from zero, indicating through Eq. [3] an elongation of the potential. The 
weighted mean elongation is (e pot sin(2<yS2)} = -0.067 ± 0.085, which is the la rgest elongation we h ave measured in our sample. 
It is due to its large uncertainty consistent with both the CDM predictions by [Havashi et alj d2007l) and with a round potential. 
The large uncertainty for the elongation arises from the relatively large error bars of the inclination values, which enter into the 
uncertainty in e pot as a fourth power (cf. Eq.0. The median amplitude of each harmonic component is fairly low (A m < 4 kms" 1 ), 
both for the entire galaxy and for the inner 1 kpc. The distribution of A r (r) shows that the amplitude of non-circular motions is 
small, especially in the inner parts of NGC 7793. Its median value (averaged over the entire galaxy) is A r ~ 5 km s" 1 . The median 
of the absolute residual velocity field is 2.2 kms" 1 , again showing that a harmonic expansion up to third order did cover most of 
the non-circular motions. 

B. THE ATLAS 

The center estimates and the harmonic decompositions are shown in Figs. 12511421 Following is a description of the content and 
layout of the figures: Each figure consists of four rows. From top to bottom, they contain: 

Row 1: The radial variation of the center position. Shown is the offset from the pointing center in arcsecond. The top panel 
shows the offsets in the X (or right ascension) direction, with positive 8X to the west; the bottom panel those in the Y (declination) 
direction, with positive SY to the north. The filled circles represent the individual center positions derived using ROTCUR, and 
the error bars indicate the formal uncertainties derived by ROTCUR. The two dotted, vertical lines indicate the radial range over 
which the positions were averaged in order to derive a kinematic center. The kinematic center is indicated by the solid horizontal 
line, and the la standard deviation by the dotted, horizontal lines. For those cases where the kinematic center was not used as 
our best center position, we additionally show the best estimate with a dashed, horizontal line. 

Row 2: The inner 150" x 150" of the galaxy. The total intensity H I map is shown in grayscales, and the beam size is indicated 
in the bottom-left corner. The black contours a re drawn from the velocity field which was used for th e estimate of the dynamical 
center (bulk velocity field from lOh et alj2008l for IC 2574, h ermite velocity field from lde Blok et al. 2008 for the other galaxies). 
The thick line represents the systemic velocity as derived in de Blo k et all (120081) . The thin black contours overlaid on the thick 
white contours belong to the 3.6 /im IRAC image and are usually given at a 2, 5, 10, 20, and 50 percent level of the maximum 
flux in the image. The white contours are taken from the THINGS radio continuum maps and are usually given at a 10, 20, and 
50 percent level of the maximum flux. The black, filled circles indicate the individual center positions from ROTCUR and the 
black cross represents the chosen dynamical center and its uncertainty. The derived center from the 3.6 /im image is shown as a 
gray, filled triangle, and the one from the radio continuum as a black, open triangle. 

Inset: To better highlight the different center estimates, we additionally show in the upper-right corner an inset with only the 
central few arcseconds. For clarity, we omit the total H I map and the velocity field contours in the inset. The contours from the 
3.6 /im image are shown in black and are given at the same intensity levels as in the main plot. The same holds for the radio 
continuum contours, which are shown here in gray. The individual center estimates from ROTCUR are shown as small crosses. In 
the inset, the beam size is indicated by the thick dashed ellipse, which is centered on our best center position. 

Row 3: The error bars of all plots in this row are the formal uncertainties as derived by RESWRI. 
Left panel: This panel consists of six sub-panels, all of which are plotted against radius in arcsecond. The upper-left sub-panel 
shows ci, the amplitude of the circular velocity. The other five sub-panels show the amplitudes of the non-circular components, 
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namely the second and third order component of the cosine term (02,03) and the first, second, and third order component of the 
sine term (si,S2,Ss). The axis scale of these five sub-panels is generally -20 to +20 kms" 1 . The solid, vertical line in the panel 
containing c\ indicates the radius corresponding to 1 kpc. All quantities shown in Row 3 are corrected for inclination effects. 
Right panel: From top to bottom: systemic velocity cq; inclination angle z; and position angle PA, all plotted against radius 
in arcsecond. The dashed, horizontal line in all three sub-panels indicates the weighted mean, using the inverse square of the 
uncertainties as weight. 

Row 4: Left panel: This panel consists of two sub-panels. The left sub-panel shows A m , the median amplitudes for each 
harmonic order m. The values are calculated by radially averaging the quadratically added amplitudes as described by Eqs. [4] 
and [5] We show the median as determined for the entire extent of the galaxy (filled circles) as well as using only the inner 1 
kpc (open circles). The error bars indicate the lower and upper quartile respectively. The right sub-panel shows A r (r), the radial 
distribution of the amplitude of all non-circular motions derived following Eq.|6] The error bars of A r (r) are derived with formal 
error propagation assuming a Gaussian error distribution. Both A m and A,(r) are corrected for inclination effects. 
Right panel: This panel contains the elongation of the potential (e pot sin(2(^2)} (cf., Eq. [3]) vs. radius. The elongation is also 
corrected for inclination. The solid, black line indicates the zero level, the dotted gray one represents the weighted mean elon- 
gation, and the dashed, gray ones its standard deviation. The error bars shown here are derived using a formal error propagation 
assuming a Gaussian error distribution. 
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FIG. 25. — Summary panel for NGC 925. Lines and symbols are described in the t ext, A ppendix IBl The IRAC center shown was derived using ELLFIT. We 
omitted the IRAC contours at the 2% and 5% levels for clarity reasons. See Appendix |A. 1 | for a discussion of this galaxy. 
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FIG. 26. — Summary panel for NGC 2366. Lines and symbols are described in the text, Appendix IBl The IRAC center shown was derived using ELLFIT, and 
the contours of the IRAC image are given at the 20% and 50% level of the maximum intensity. See Appendix |A.2| for a discussion of this galaxy. 
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FIG. 27. — Summary panel for NGC 2 403. L ines and symbols are described in the text, Appendix IBl The IRAC contours are given at the 20, 50, and 80% level 
of the maximum intensity. See Appendix |A.3| for a discussion of this galaxy. 
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FIG. 28. — Summary panel for NGC 2841. Lines and symbols are described in the text, Appendix|B] See Appendix |A.4| for a discussion of this galaxy. 
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FIG. 29. — Summary panel for NGC 2903. Lines and symbols are described in the text, Appendix|B] See Appendix IA.5 I for a discussion of this galaxy. 
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FIG. 30. — Summary panel for NGC 2976. Lines and symbols are described in the text, Appendix[B] The IRAC contours are given at the 10, 20, 50, and 80% 
level of the maximum intensity. See Appendix |A.6| for a discussion of this galaxy. 
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FIG. 3 1 . — Summary panel for NGC 3031. Lines and symbols are described in the text, Appendix|B] The bright radio continuum source which is s till vi sible 
in the total intensity H I map was blanked out for clarity reasons. Due to the central H I deficiency, no velocity contours are visible. See Appendix |A.7| for a 
discussion of this galaxy. 
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FIG. 32. — Summary panel for IC 2574. Lines and symbols are described in the tex t, Ap pendix IBl No estimate could be derived from the IRAC and radio 
continuum images, and their respective contours are therefore not shown. See Appendix |A.9| for a discussion of this galaxy. 
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FIG. 33 . — Summary panel for NGC 3521 . Lines and symbols are described in the text, Appendix[B] See Appendix IA. 1 Ol for a discussion of this galaxy. 
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FIG. 34. — Summary panel for NGC 3621. Lines and symbols are described in the text, Appendix IBl The IRAC contours are given at the 20, 50, and 80% level 
of the maximum intensity. Only data with r < 600" was used for radial averaging of A,„ or e p oi. See Appendix lA. 1 1 I for a discussion of this galaxy. 
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FIG. 35. — Summary panel for NGC 3627. Lines and symbols are described in the text, Appendix|B] The axis scale on the panels showing the non-circular 
components runs from -30 to 40 kms~' . See Appendix lA. 12l for a discussion of this galaxy. 
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FIG. 36. — Summary panel for NGC 4736. Lines and symbols are described in the text, Appendix [B] The IRAC contours are given at the 2, 10, 50, and 
80% level of the maximum intensity. The axis scale on the panels showing the non-circular components runs from —30 to 35 kms~'. See Appendix IA. 1 3 I for a 
discussion of this galaxy. 
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FIG. 37. — Summary panel for DDO 154. Lines and symbols are described in the te xt, App endix 151 No estimate could be derived from the IRAC and radio 
continuum images, and their respective contours are therefore not shown. See Appendix I A. 1 4 [ for a discussion of this galaxy. 
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FIG. 38. — Summary panel for NGC 4826. Lines and symbols are described in t he text , Appendix 151 The lower panel is missing given that we were not able 
to derive a meaningful harmonic decomposition of the velocity field. See Appendix lA. 1 5l for a discussion of this galaxy. 
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FIG. 39. — Summary panel for NGC 5055. Lines and symbols are described in the text, Appendix|B] The axis scale on the panels showing the non-circular 
components runs from -30 to 30 kms~' . Only data with r < 450" was used for radial averaging of A m or e pot . See Appendix l A. 1 6l for a discussion of this galaxy. 
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FIG. 40. — Summary panel for NGC 6946. Lines and symbols are described in the text, Appendix[B] The IRAC contours are given at the 2, 5, and 50% level 
of the maximum intensity. The axis scale on the panels showing the non-circular components runs from —45 to 25 kms~' . Only data with r < 420" was used for 
radial averaging of A„, or e pot . See Appendix lA. 17l for a discussion of this galaxy. 
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Fig. 41.— 



Summary panel for NGC 7331. Lines and symbols are described in the text, Appendix[B] See Appendix IA. 1 8l for a discussion of this galaxy. 
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FIG. 42. — Summ ary panel for NGC 7793. Lines and symbols are described in the text, Appendix IBl The 2% IRAC contour was omitted for clarity reasons. 
See Appendix |A. 19| for a discussion of this galaxy. 



